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FOREWORD 


The  investigation  reported  herein  is  part  of  the  research 
project  at  The  Ohio  State  University,  Columbus,  Ohio  sup¬ 
ported  by  the  Air  Force  Office  of  Scientific  Research  Grant 
83-00-55.  Lt.  Col.  John  J.  Allen  was  the  Program  Manager  at 
the  commencement  of  the  Project.  Lt.  Col.  Lawrence  D.  Hok- 
anson  was  the  Program  Manager  from  July  1,1983.  The  Research 
project  was  started  in  February  1,  1983  and  is  continuing. 
The  present  report  documents  part  of  the  work  done  up  to 
January  31,  1984.  At  The  Ohio  State  University,  the  Project 
is  being  supervised  by  Or.  Ranbir  S.  Sandhu,  Professor,  De¬ 
partment  of  Civil  Engineering.  The  computer  program  develop¬ 
ment  and  the  analyses  reported  herein  were  carried  out  by 
Mahantesh  S.  Hiremath,  Graduate  Research  Associate.  The  In¬ 
struction  and  Research  Computer  Center  at  The  Ohio  State 
University  provided  the  computational  facilities. 


ABSTRACT 


The  'Engineering  Approach'  to  dynamic  response  ana¬ 
lysis  of  saturated  sand  deposits  is  investigated.  The  logic 
on  which  the  methodology  is  based  is  described  and  imple¬ 
mented  in  a  computer  program.  Application  of  the  procedure 
to  obtain  dynamic  response  of  a  saturated  sand  layer,  in¬ 
cluding  pore  pressure,  shear  stress  and  acceleration  varia¬ 
tions  under  two  different  ground  excitation  histories  is 
given.  Two  alternative  numerical  procedures  are  investigat¬ 
ed.  The  results  are  compared  with  those  reported  by  Finn. 
Limitations  of  the  approach  are  discussed. 
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SECTION  I 


INTRODUCTION 

Since  the  early  forties,  the  mechanism  underyling  li¬ 
quefaction  of  saturated  sands  subjected  to  dynamic  loading 
has  attracted  the  attention  of  engineers  and  geologists. 
Particularly,  the  events  following  the  earthquakes  at  Alaska 
and  Niigata  in  1964  and  the  extensive  damage  caused  by  li¬ 
quefaction  of  soil  in  both  the  cases  led  to  extensive  geo¬ 
technical  engineering  effort  to  study  of  response  of  satu¬ 
rated  sand  deposits  to  earthquake  motions.  Along  with  this, 
the  need  to  consider  this  problem  in  design  of  important 
structures  such  as  dams,  nuclear  reactors  etc.  has  played  a 
major  role  in  continued  interest  in  this  area. 

The  effect  of  earthquake  induced  shear  stresses  and 
shear  strains  is  to  cause  slip  at  grain  to  grain  contact, 
which  results  in  volumetric  compaction  in  dry  sands.  "In  sa¬ 
turated  sands  however,  any  volume  change  of  sand  implies  a 
corresponding  quantity  of  fluid  flowing  out.  This  is  neces¬ 
sarily  a  time  depedent  process.  Hence  volumetric  compaction 
takes  place  gradually  by  drainage.  According  to  Finn  [6], 
the  mechanism  of  deformation  and  pore  water  pressure  build¬ 
up  is  as  follows:  In  saturated  sands, the  volumetric  compac- 


tion  is  retarded  because  the  water  cannot  drain  instantane¬ 
ously  to  accomodate  the  volume  change.  Consequently,  the  re¬ 
laxing  sand  skeleton  transfers  some  of  its  i ntergranul ar  or 
effective  stresses  to  the  pore  water  and  the  pore  water 
pressure  increases.  The  corresponding  reduction  in  effective 
stress  leads  to  a  structural  rebound  in  the  sand  skeleton  to 
absorb  the  difference  in  volume  between  the  compaction  due 
to  the  grain  slips  and  the  reduction  in  pore  water  volume 
due  to  increased  pore  water  pressure  and  drainage.  In  the 
extreme  case,  the  pore  water  pressure  developed  during  the 
earthquake  may  increase  until  all  the  i ntergranul ar  or  ef¬ 
fective  stress  has  been  eliminated.  In  this  state,  the  sand 
has  no  significant  shearing  resistance  and  deforms  like  a 
liquid.  Liquefaction  is  said  to  have  occured". 

The  early  attempts  at  determining  the  response  to 
cyclic  loading  of  a  horizontally  layered  sand  deposit  were 
based  on  total  stress  procedures.  The  volume  changes  and  the 
porewater  pressures  that  develop  during  cyclic  loading  were 
assumed  to  depend  upon  the  shear  strains.  The  dynamic  shear 
strains  in  turn  depended  upon  the  stiffness  and  damping 
characteristics  of  the  sand  layers.  For  a  horizontal  shear 
wave  propagating  vertically  through  a  horizontally  layered 
saturated  sand  deposit,  assuming  that  the  sand  undergoes 
simple  shear  deformations  only,  the  stiffness  at  any  time 
may  be  expressed  as  a  function  of  the  shear  modulus.  In  fun- 


damental  studies  of  effective  stress-strain  relations  for 
sands,  Seed  and  Idriss  [1]  and  Hardin  and  D  me  v  i  ch  [2] 
showed  the  shear  modulus  to  be  function  of  the  mean  effec¬ 
tive  stress  and  the  shear  strain.  The  initial  effective 
stresses  were  used  in  computation  of  the  initial  value  of 
the  shear  modulus.  Subsequently,  the  dependence  on  shear 
strains  was  allowed  for.  However,  increase  of  the  porewater 
pressure  during  shaking  and  consequent  decrease  in  the  ef¬ 
fective  stresses  in  the  layer,  which  in  turn  decrease  the 
mean  normal  effective  stresses,  was  not  taken  into  account 
in  calculating  the  shear  modulus. 

An  improved  method  for  problems  of  liquefaction  in  two 
dimensional  systems  using  iterative  approach  was  presented 
by  Seed,  Lee  and  Idrtss  [3].  In  study  of  Sheffield  dam,  they 
assigned  zero  shear  modulus  to  those  regions  of  the  dam  that 
liquefied.  Proceeding  in  this  manner,  it  was  possible  to 
predict  the  progressive  developement  of  liquefaction.  Howev¬ 
er,  in  case  of  one-dimensional  problems,  since  liquefaction 
is  supposed  to  occur  at  the  same  time  at  points  on  the  same 
level,  the  above  method  could  not  be  applied.  Also,  effect 
of  the  porewater  pressure  on  the  shear  modulus  prior  to  li¬ 
quefaction,  was  still  not  accounted  for. 

Streeter  et  al  [4]  illustrated  a  numerical  method  of 
dynamic  response,  that  accounts  for  the  effect  of  pore  pres- 


sure.  The  response  of  the  pore  water  and  the  granular  skele¬ 
ton  were  treated  separately  as  uncoupled  problems.  Pore 
pressure  effects  were  incorporated  by  specifying  pseudocycl- 
ic  and  permanent  volume  changes.  However,  a  method  for  as¬ 
signing  values  to  such  volume  changes  was  not  established. 

Liou,  Streeter  and  Richart  [5]  introduced  a  different 
model  for  the  development  of  pore  pressure.  A  fixed  rela¬ 
tionship  between  shear  modulus  and  the  constrained  rebound 
modulus  was  assumed  to  exist.  The  adequacy  of  this  model  is 
yet  to  be  established. 

It  would  seem  reasonable  to  expect  that  a  realistic 
representati on  of  the  dynamic  response  of  saturated  sand  de¬ 
posit  would  require  allowing  for  progressive  development  of 
pore  pressure  simultaneously  with  the  dynamic  analysis  of 
displacements.  A  method  of  effective  stress  analysis  was 
presented  by  Finn  et  al  [6].  In  this,  variation  in  pore  wa¬ 
ter  pressure  and  strains  was  allowed  by  periodic  updates. 
The  constitutive  laws  account  for  important  factors  known  to 
affect  the  response  characteristic  of  saturated  sands  to  dy¬ 
namic  loading,  including  simultaneous  generation  and  dissi¬ 
pation  of  pore  pressure.  This  method  and  its  variants  are 
the  current  state-of-the-art.  As  part  of  the  present  re¬ 
search  effort,  this  'engineering  approach'  was  implemented 
in  a  computer  program.  Specific  application  considered  was 


horizontally  layered  saturated  sand  deposits  subjected  to  a 
horizontal  shear  wave  applied  at  the  base  and  propagating 
vertically  up.  This  theory  does  not  allow  for  continuous 
and  simultaneous  changes  in  pore  water  pressure  and  strain. 
Periodic  updates  are  an  improvement  over  no  updates  at  all 
but  as  Christian  [24]  noted  for  the  case  of  consolidation  of 
soils,  the  sequential  solution  of  two  differential  equations 
did  not  give  correct  solutions  and  was  very  sensitive  to 
mesh  size.  For  the  case  of  dynamic  response  analysis,  the 
need  to  be  careful  is  even  greater.  The  sequence  of  occur¬ 
rences  assumed  in  the  'engineering  approach'  viz.  shear 
stresses  cause  volume  change,  volume  change  results  in  pore 
water  pressure  changes,  pore  pressure  dissipation  follows, 
pore  pressures  determine  effective  stresses,  which  in  turn 
define  shear  modulus  to  be  used  subsequently  represents  the 
physical  phenomenon  only  approximately.  The  theory  has  been 
attractive  because  of  its  'simplicity'  and  the  fact  that  no 
alternatives  were  available.  We  note  that  the  theory  is  far 
from  simple  insomuch  as  the  procedure  depends  upon  insitu 
condition  of  soil  deposit  and  uses  empirical  relations  based 
on  experimentally  obtained  soil  properties.  The  procedure  is 
site-specific.  The  numerical  solution  scheme  is  not  quite 
reliable.  The  results  of  the  analysis  can  thus  be  only  qua¬ 
litative.  The  approach  cannot  determine  the  post-liquefac¬ 
tion  pore  pressures  correctly.  It  has  been  used  only  for 
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one-dimensional  wave  propagation  and  cannot  be  extended  to 
two  and  t h ree-d i mens i ona 1  cases.  Use  of  this  theory  requires 
considerable  experience  and  "judgement''  in  addition  to  ex¬ 
tensive  laboratory  testing  program  to  get  useful  results. 


Biot's  [20]  theory  of  coupled  motion  and  fluid  flow 
has  been  used  by  Ghaboussi  and  Wilson  [23]  to  develop  compu¬ 
ter  program  for  dynamic  response  analysis  of  soils.  Howev¬ 
er,  the  procedures  are  applicable  only  to  linear  elastic 
soils,  use  several  other  assumptions  and  have  not  been  veri¬ 
fied  against  either  theoretical  solutions  or  experimental 
data  in  the  laboratory  or  in  the  field.  An  investigation  of 
that  approach  is  the  subject  of  a  separate  report. 


SECTION  II  covers  the  theoretical  background  for  the 
'engineering  approach1  including  constitutive  laws,  direct 
nonlinear  analysis,  coupled  dynamic  analysis  and  pore  pres¬ 
sure  response.  SECTION  III  describes  the  example  problem. 
Results  of  the  analysis  are  given  in  SECTION  IV.  SECTION  V 
contains  the  conclusions  of  the  present  study.  The  computer 
program  description,  instruction  for  usage  and  Fortran  list¬ 
ing  along  with  sample  input  and  output  are  given  in  the  Ap- 
pendi ces . 


SECTION  II 


THEORETICAL  FORMULATION 

2.1  PHYSICAL  INTERPRETATION  FOR  DYNAMIC  RESPONSE 

It  is  generally  recognized  that  the  basic  cause  of  li¬ 
quefaction  of  saturated  cohesionless  soils  subjected  to  dy¬ 
namic  loading  is  the  build-up  of  excess  hydrostatic  pres¬ 
sure.  These  stresses  could  be  due  to  upward  propagation  of 
shear  waves  associated  with  a  seismic  event  or  due  to  blast 
or  explosion. 

In  the  'engineering  approach',  the  basic  line  of  rea¬ 
soning  is  as  follows.  As  a  consequence  of  dynamic  stresses, 
caused  by  applied  loading,  the  sand  tends  to  become  more 
compact  with  transfer  of  stress  to  the  pore  water.  As  a  re¬ 
sult,  the  sand  grain  structure  rebounds  to  the  extent  re¬ 
quired  to  keep  the  volume  constant  and  this  interplay  of  vo¬ 
lume  reduction  and  soil  structure  rebound  determines  the 
magnitude  of  the  increase  in  porewater  pressure  in  the  sand 
deposit.  Under  reduced  effective  stress  and  continued  motion 
slips  at  the  grain  contacts  are  possible  resulting  in  a  more 
stable  sand  structure  under  the  reduced  effective  stress 


[7].  Actually,  the  changes  in  effective  stress,  skeletal 
structure  and  pore  water  pressure  occur  s i mu  1 taneous ly  rath¬ 
er  than  as  a  sequence  of  events  listed  above.  In  a  more  com¬ 
plete  theory,  these  occurrences  would  be  treated  as  simulta¬ 
neous  . 

In  the  'engineering  approach',  to  compute  the  response 
of  saturated  sand  layers  under  applied  shear  excitation,  the 
following  factors  need  to  be  considered: 

1.  Initial  shear  modulus  in  situ. 

2.  Variation  of  shear  modulus  with  shear  strain. 

3.  Generation  of  pore  water  pressure  due  to  straining  of 
soil  and  dissipation  thereof  due  to  spatial  flow. 

4.  Changes  in  effective  mean  stress. 

5.  Damping  and  its  influence.. 

6.  Processes  leading  to  increased  resistance  termed  col¬ 
lectively  as  hardening. 

2.2  CONSTITUTIVE  RELATIONS 

Constitutive  relations  for  dry  or  saturated  sands  in 
simple  shear  were  presented  by  Finn[8].  These  covered  the 
factors  cited  earlier  and  subsequent  analytical  procedures 
included  analysis  of  dynamic  response  of  saturated  sands  un¬ 
der  field  loading  conditions  along  with  the  dissipation  of 
pore  pressure  [’9,10]. 


The  initial  response  of  sand  to  dynamic  loading  depen¬ 
ds  upon  its  state  in  situ,  as  indicated  by  its  shear  modulus 

Finn  [10]  assumed  that  upto  the  point  of  first  reversal 
in  loading,  the  response  of  the  sand  follows  the  hyperbolic 
stress-strain  relationship  formulated  by  Kondner  and  Zelasko 
[11]  and  also  adopted  by  Hardin  and  Drnevich  [2].  This  re¬ 
lationship  is  represented  by  the  following  equation: 

T.«  f(  Y  ) 

or  r  .  Gn0v/  (l  ♦  Gm0  y,  Tm0)  (1) 

in  which  r  is  the  shear  stress  at  the  strain  amplitude  of  Y  , 
Gm0  is  the  initial  maximum  tangent  modulus  and  rmQ  is  the 
maximum  shear  that  can  be  applied  to  the  sand  in  its  initial 
state  without  failure  as  shown  in  Figure  (1). 

The  quantities  G^g  and  rmg  were  given  in  units  of  po¬ 
unds  per  square  foot  by  equations  proposed  by  Hardin  [7], 
These  were  modified  for  horizontal  sand  layer  as  [6] 

Gm0  =  14760  [(2.973-e)2/(l+e)]  [ ( l  +  2k Q ) / 2  ] 1 ' 2  -[<  \  (2) 

Tm0  =  [[0.5(  l+k0)  sin  *]2  -  [0.5(l-kg)]2]1/2  <'v  (3) 

in  which  e  =  void  ratio  a-nd  is  limited  to  values  less  than 
2,  <'  =  vertical  effective  stress  in  psf,  Kn  =  coefficient 


of  earth  pressure  at  rest  and  ^  =  effective  angle  of  shear 
resistance.  For  unloading  or  reloading  phase,  the  stress- 
strain  curve  is  described  by  Masing's  [21]  rule  i.  e. 

(  T  -  rp)/2  =  f c (  y-  y r)/2]  (4) 

where  {  r  ,  y  )  is  the  point  at  which  loading  reversal  occu¬ 
rs.  Equation  (1)  through  Equation  (4)  imply  that  the  she¬ 
ar  modulus  is  strain-dependent  and  damping  is  hysteretic  [8] 

2.3  HARDENING  AND  PORE  PRESSURE  EFFECTS 

During  cyclic  loading  slippage  at  the  grain  contacts 
results  in  volumetric  compaction  [7,12,13]  and  increased  va¬ 
lue  of  kg,  the  coefficient  of  effective  lateral  stress  [14]. 
Both  these  effects  stiffen  the  sand  against  further  deforma¬ 
tion.  It  is  probable  that  slips  at  the  grain  contacts  result 
in  a  more  stable  structure  under  given  effective  stress  [7]. 
These  processes  leading  to  increased  resistance  are  referred 
to  collectively  as  hardening.  Martin  [12]  suggested 
a  stress-strain  relationship  in  the  form 

Thv  a  yfc'v  /{  a  +  by  }  (5) 

in  which  Thv  and  y  are  associated  horizontal  shear  stress 
and  shear  strain  respectively  and  <'  is  the  current  effect- 


i ve  stress . 


The  expression  for  maximum  shear  modulus  G  amd  maxi¬ 
mum  shear  stress  rmn  in  terms  of  the  initial  values  Gmg  and 
t^q  and  volumetric  strain  are  given  as 

G,n  *  Gm0  [  1  *  «,d  '  I  H1  *  h2  «vd  >  J  (6) 

Tmn  ■  Tm0  C  1  *  «vd  '  <  H3  *  H4  «,d  >  J  <7> 

where  6yd  is  the  accumulated  volumetric  strain  and  H^,  H2, 

H3  and  are  constants. 

2.4  GENERATION  OF  PORE  PRESSURE 

According  to  Martin  [12],  during  an  undrained  shear 
test,  the  cyclic  shear  strain  causes  an  increase  in  pore  wa¬ 
ter  pressure,  Au,  expressed  as 

An  ’At  v<j  /[  1/E  r  *  np/Kw  ]  (8) 

in  which  Er  *  one-dimensional  rebound  modulus  of  sand  at  an 
effective  stress  4y  and  Kw  is  the  bulk  modulus  of  water.  For 
saturated  sands  K  >>  E  and  therefore  Equation  (8)  becomes 


Finn  [8]  expected  A«vd  to  be  a  function  of  total  accumulated 
volumetric  strain,  and  amplitude  of  the  shear  strain  and 
proposed  the  relationship 


-  c,(  y-  c2  €vd)  +  c,  €uA/(  y  +  c,€uh) 


3  vd' 


'4  vd 


For  crystal  silica  sand  with  relative  density  Dr  =  45^,  Kon- 
dner  [11]  gave  the  values  of  the  constants  as  :  C  ^  =  0.80, 

C2  ■  0.79,  C3  =  0.45  and  C4  =  0.63,  when  €vd  and  are 
expressed  as  percentages.  An  expression  for  the  rebound  mod¬ 
ulus*  at  any  effective  stress  level,  <*v ,  was  given  by 
Martin  [12]  as 


V  *  (  <*  )*•■  /  m  K,  (  <'  )n-m 


in  which  ^ wq  =  initial  value  of  the  effective  stress  and 
K2,  m  and  n  are  experimental  constants.  For  crystal  silica 
sand  at  =  45%,  the  values  of  K2»  m  and  n  were  given  by 
Kondner  [11]  as  0.0025,  0.43  and  0.62  respectively. 

The  increment  in  pore  pressure,  Au,  during  a  given 
loading  cycle  with  maximum  shear  strain,  y,  may  now  be  com¬ 
puted  using  Equation  (9)  through  Equation  (11).  For  satu¬ 
rated  sands,  the  maximum  shear  moduli  and  maximum  allowable 
shear  stresses  for  nth  loading  cycle  are  related  to  initial 
val ues  by  [6] 


Gmn  =  GmO^  +  €»n/(  Hi  +  <ufi> 


vd'  1  nl  T  n2  c  vd  '  -^v'  N  vO 


^mn  ~TmO^  +  €vd^  ^*3  +  ^4  *vd^^v/l 


in  which  <^q  is  the  effective  vertical  stress  and  <*v  is  the 
effective  vertical  stress  at  the  beginning  of  the  n  th  cycle 


2.5  DISSIPATION  OF  PORE  WATER  PRESSURE 


If  the  saturated  sand  layer  can  drain  during  motion, 
there  will  be  simultaneous  generation  and  dissipation  of 
porewater  pressure.  Thus,  the  rate  of  increase  of  porewater 
pressure  will  be  less  than  for  completely  undrained  sand. 
The  distribution  of  porewater  pressure  at  time  is  given  by 
the  equation 


U,t  ”  ^K/,yw^u,z^,z  +  ^r  ^  €vd-^,t 


(14) 


in  which  a  subscripted  comma  indicates  the  derivative  with 
respect  to  the  variable(s)  following  it,  u  is  the  pore  wat¬ 
er  pressure,  K  is  the  coefficient  of  permeability  and  y is 
the  unit  weight  of  water.  The  term  containing  represents 
the  internal  generation  of  pore  water  pressure  and  is  simi¬ 
lar  to  a  source  term  in  a  heat  conduction  equation. 


The  porewater  pressure  any  time  is  in  accordance  with 
Equation  (14)  and  reflects  the  net  effects  of  contemporane¬ 
ous  generation  due  to  volume  change  and  redistribution  due 
to  drainage. 

The  shear  modulus,  G,  expressed  in  terms  of  mean  nor¬ 
mal  effective  stress  ^  was  given  by  Seed  [15]  as 

G  =  1000  K2  (  <'m)1/2  (15) 

in  which  <2  is  a  parameter  that  varies  with  shear  strain  as 
shown  in  Figure  (2)  for  relative  density  =  0.45. 

Equation  (14)  must  be  solved  in  conjunction  with  the 
equation  of  motion  of  the  sand  layer  in  order  to  continually 
update  the  values  of  the  p-ore  water  pressure  during  dynamic 
loading. 

2.6  COUPLED  DYNAMIC  ANALYSIS 

Consider  a  stratum  of  saturated  sand  of  infinite  la¬ 
teral  extent  resting  on  horizontal  bedrock.  The  properties 
of  the  stratum  may  vary  only  in  vertical  direction.  Let  the 
applied  excitation  consist  of  horizontal  shear  waves  applied 
at  the  base  propagating  vertically  upwards.  The  dynamic  res¬ 
ponse  in  this  case,  in  the  'engineering  approach'  is  be  con¬ 
sidered  equivalent  to  that  of  a  shear  beam. 


v.v 
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2.6.1  Differential  Equation  of  Motion 


Since  the  properties  of  the  soil  deposit  may  vary  only 


along  the  vertical,  the  horizontal  sand  deposit  is  divided 


into  layers  as  shown  in  Figure  (3). 


The  layered  system  is  modelled  as  a  lumped  mass  system 


with  lumped  masses  (m^  m  .  mn).  The 


mass 


Yx  hj  /  (2g) 


(16) 


is  assigned  to  the  top  of  the  soil  layer  No.l.  Here  and 


are  the  unit  weight  and  the  thickness,  respect i ve 1 y ,  of  soil 


layer  No.  1.  The  other  lumped  masses  associated  with  inter¬ 


layer  locations  are 


'l  *  C  >Vl  hf-1  *  yi  "i  3/(2 9) 


(H) 


where  i  -  1,2, 


The  masses  are  assumed  to  be  connected  by  non-linear 


springs  having  initial  properties  given  by  Equation  (1)  and 


modified  by  using  Equations  (14)  and  (15).  These  equations 


reflect  non-linear,  strain-dependent,  hysteretic  behaviour 


of  sand.  The  spring  constants  are  given  by 


Ki  ■  v/  "i 


*<  •*.  *  .  .  -*  ■*. 


*1  =  Gi  '  h1 


(18) 


is  the  spring  constant  of  the  spring  connecting  masses 
and  and  is  the  shear  modulus  of  the  layer  i.  The 

equation  of  motion  of  the  system  can  be  writen  as  [16], 


CM]  {u>  ♦  [C]  0}  +  CK]  =  {R(t)> 


(19) 


['MJ  is  the  diagonal  mass  matrix,  [C]  is  the  viscous  damping 
matrix,  [K]  is  the  tridiagonal  symmetric  non-linear  stiff¬ 
ness  matrix,  {R(t)>  is  the  excitation  applied  to  the  system 
and  {u},  {u}  and  {u}  are  the  displacements,  velocities 
and  accelerations  respectively  of  the  masses.  Thus, 


diag  CMJ  =  (  mp  m2> 


•  %> 


(20) 


and  [K]  is  a  tridiagonal  matrix  with 


K11  *  K1 


Ki_1  +  Ki  for  i * j 


K  .  .  s  -  K 

IT  'M 


for  i = j - 1 


K  .  .  =  -  K  . 

1  J  J  1 


for  i =j  +1 


K(j  .  0. 


for  al 1  other  i  ,  j . 


The  load  vector,  for  the  case  of  acceleration  applied 

to  the  base  of  the  system,  is  given  by 


{R(t)>  =  -  col  (mlf  m2. 


Cug> 


To  set  up  the  damping  matrix,  which  would  allow  uncoupling 
of  the  damping  forces  in  the  same  way  as  the  uncoupling  of 
the  inertia  and  elastic  forces,  Caughey  [25]  proposed  use  of 
a  finite  series  in  following  form 


C  »  M  a  [m'1  K]b 


(23) 


where  b  is  number  of  terms  in  the  series  and  a  is  proportio¬ 
nality  factor.  For  b  =  1,  this  reduces  to  Rayleigh  damping 
with  C  being  a  linear  combination  of  M  and  K.  In  general,  d- 
amping  ratio  associated  with  any  mode  n  is  given  by  [16] 


fn  =  ‘/(2  »n>  ah 


(24) 


A  special  case  of  Rayleigh  damping  is  the  stiffness  proporti 
onal  damping  in  which 

[C]  *  a :  [ « ]  (25) 

Here  the  higher  modes  of  the  structure  are  very  heavily 
damped,  since  damping  ratio  is  directly  proportional  to  fre¬ 
quency.  The  fundamental  natural  frequency  of  a  cohension- 
less  soil  strata  of  depth  H,  shear  modulus  G  and  unit 
weight  of  X  is  given  by  [17] 

wn  *  1.7510  «/g /( 1.2)  (H)S/6  (26) 

In  the  present  analysis,  the  natural  frequency  of  the  stra¬ 
tum  was  calculated  taking  the  values  of  H,  G  and  of  the  sa¬ 
turated  portion  only  using  Equation  (26)  and  then  it  was 
used  in  Equation  (24)  to  calculate  a^. 

2.7  SOLUTION  OF  THE  DYNAMIC  EQUATIONS  OF  MOTION 

Following  Finn  [6],  Equation  (19)  was  solved  numeri¬ 
cally  using  Newmark's  unconditionally  stable  integration  op¬ 
erator  [18].  For  completeness,  the  method  is  described 
briefly  in  the  following  paragraphs. 


2.7.1  Newmark's  Method 


If  u,  u  and  u  denote  the  displacements,  velocities  and 
accelerations  respectively,  then  Newmark's  method  assumes. 


u  ( t  +At )  *  u  ( t )  +  [  ( 1  -  8  )  u  ( t )  +  5u(t+At)]At 


u  ( t  +  A  t )  s  u(t)  +  ii  ( t )  A  t  +  [  (  0 . 5  -  a  )  u  ( t )  +  «u(t+At)At 


(28) 


where  S  and  a  are  parameters  that  can  be  determined  to  ob¬ 
tain  integration  accuracy  and  stability.  Newmark  [18]  pro¬ 
posed  an  unconditionally  stable  scheme  of  constant-average- 
acceleration,  in  which  case  5  =  0.50  and  a-  0.25.  The 
equilibrium  Equation  (19)  at  time  (t+At)  is  used  to  evalu¬ 
ate  u(t+ At)  i  -e. 


[M]  u  ( t  +  At)  +  [C]  u(t+  At)  +  [K]  u  ( t+  A  t )  «  (R(ttAt)l  (29) 


From  Equation  (28),  u(t  +  At)  is, 


u(t+At)  =  [l/(«At2)]  [u(t+ At)  -  u(t)  -  u  ( t )  At] 


(0.5a  -  1)  u(t) 


23 

Substituting  this  in  Equation  (27),  an  expression  for  u(t+ 
At)  is  obtained  in  terms  of  u(t+At)  as>  • 

u(t+  At)  *  [1  -  5 /a  ]  u(t)  +  [(2  -  5/a  )  At/2]  u(t) 

# 

+  [  5 /a  At]  [u(t+  At)  -  u ( t )  ] 


These  two  equations  can  be  used  in  Equation  (29)  to  >;t 
f ol lowing 


u  ( t  +  A  t )  [  K  +  M  (1/aAt2)  +  C  (  5  /a  At )  ] 


*  R(  t+  At )  +  M  [  ( 1  /  a  At 2  )  u  ( t )  +  (1/a  At)u(t)  +  ( 1/2 «  -1  )u  (t )] 

+  C  [(  5  /  aAt )  u  ( t )  +  (5/a  -  1)  G(t)  +  (  At/2)  (5/a  -  2)  u(t)] 

(30) 


2.7.2  Algorithm  for  Newmark's  Method 

1.  Form  matrices  [M],  [C]  and  [K]. 

2.  Initialize  u(0),  u  ( 0 )  and  u  ( 0  ) . 

3.  Select  time  step  At,  parameters  5  and  a  such  that 
stability  conditions  [22] 

5  >  0.5 


a  *  0.25  (0.50  +  5  )2 


are  sati sf ied. 


4.  Calculate  the  coefficients 

a0  =  l/(«  At)2  *  S/(a  At  ) 

a2  =  l/  (aAt)  aj  !  0.5*1  -1 

a4  *  8/  a  -  l  a5  *  0.5  (At )  (8/  a  -  2) 

a6  -  At  (1-  5)  a  7  =  «  At 

5.  Form  effective  [K]  *  [k]  +  a g[M]  +  a^  [C] 

6.  For  each  time  step: 

a)  Calculate  effective  load  at  time  (t  +  t). 

{«"  (t+  At)}  -  {R(t+  At)3 

+  [M]  [aQ  u(t)  +  a2  u(t)  +  a  3 u  ( t )  ] 

+  [C]  [a^t)  +  a4u(t)  +  a5  u(t)] 

b)  Solve  for  displacements  at  time  (t+At). 

[£]  Cu  (t  +  At )}  *{R(t+At)} 

c)  Calculate  accelerations  and  velocities  at  time 
(t+At)  using 

u(t  +  At)  =  aQ[u(t+  At)  -  u(t)]  -  a2  u(t) 


u  ( t+  A  t )  =  u  ( t )  +  a  g  u(t)  +  a7  u(t  + At) 


(35) 


2.8  PORE  PRESSURE  DISTRIBUTION 

In  order  to  calculate  the  pore  water  pressure,  Finn 
[6]  proposed  use  of  equation  (14)  viz. 

u,t  *  Er  C(K/rw)  u>z]iZ  ♦  Er  [<vd]>t 

If  the  stratum  is  mechanically  and  hydraulically  homo¬ 
geneous,  then  K  and  f ^  are  independent  of  z.  In  that  case, 
the  above  equation  reduces  to 

\t  •  Er  ("/%>  C\223  -  Er  <36> 

Setting  D  •  Er  K/yw, 

\t  ■  0  C\22J  *  Er  C*,d3,t  <37> 

Two  classes  of  methods,  explicit  and  implicit,  are  available 
for  solution  of  Equation  (37).  In  following  sections,  we  de 
scribe  both  these  procedures. 

2.8.1  Explicit  Method 

For  a  general  diffusion  equation  of  the  type 


[•] 


Using  forward  differences  in  the  time  domain  and  central 
differences  spatially.  Equation  (38)  is  discretized  to  give 
the  explicit  relationship 


+  (D  At/h2)(.u". 


2  uj  ♦  u3.x] 


+  [  Qn  +  { 

j  +1 


-  Oj  1 


(39) 


where  the  superscripts  refer  to  the  time  step  number  and  the 
subscripts  indicate  the  number  of  the  mesh  point  in  the  spa¬ 
tial  partition.  The  difference  formula  is  illustrated  by  a 
di agram  in  Figure  (  4 ) . 


It  is  noted  that  value  of  u  at  time  Tn+1  at  point  x- 

depends  on  the  value  of  the  function  at  earlier  time  T  and 

n 

not  on  any  value  of  the  current  time.  For  stability,  this  m- 
ethod  requires  [19]  that  the  value  of  the  term  [DAt/h2]  <  .5 
where  h  denotes  the  length  of  the  spatial  interval.  Noting 
that  for  continued  saturation  by  an  incompressible  fluid  the 
flow  equals  the  volumetric  strain.  Equation  (39)  may  be  wri- 


ten  as 
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[un+1  .  un] 
J  J  J 


[D  At/h2]  [Ujn  +  1  -  2  u]  +  uj.f3 


+ 


vd  j 


n 

vd  j 


]  E, 


(40) 


or 


n+1  =  un 
J  J 

+  CD  At/h2]  [uj  +  }  -  2  +  uJ.j] 

+  Cuj  +  1  *  2  uj  +  uj-l3 


+  [ €  n+l 

vdj 


vd  j 


(41) 


Setting  [D  it/h2],  Equation  (41)  has  the  form 


»J+1-  “3  -  »$-ii 


C£r]  [A« 


(42) 


2.8.2  Implicit  Method 

The  implicit  method  using  finite  differences  called 
Crank-Nicolson  scheme  avoids  the  stability  requirement  of 
the  explicit  method  [19].  Thus,  Equation  (38)  can  be  writen 
using  central  difference  to  approximate  both  spatial  and 
temporal  derivatives  as 


-  uj  »  (DAt/2h‘)  [u'j  +  1  -  2u'!  +  uV.j 


+  “j:!  - 2  u"+1  +  uj-}] 


+  [QH  +  1 

J 


Letting  0-  (D  At/h2 ) , 


*  2(1+0 )  Ujn+i  . 

^u5  +  1  +  2(1-0  )  u]  +  i9uj-1+  AQ 


(44) 


where  Aq  denotes  (QH  +  1  -  q|? ) .  Thus,  for  finding  the  solution 

J  J 

at  time  (n+1)  requires  the  solution  of  a  set  of  linear  alge¬ 
braic  equations.  Since  the  system  is  tridiagonal,  the  cost 
per  time  step  of  this  method  is  much  greater.  But,  this 
method  is  an  order  of  magnitude  more  accurate  and  uncondi¬ 
tionally  stable.  It  permits  use  of  a  larger  time  step  re¬ 
sulting  in  significant  overall  economy,  for  comparable  accu¬ 
racy,  over  the  explicit  method.  The  computational  molecule 
of  the  method  is  shown  in  Figure  (5). 


It  is  noted  that  expression  for  over  the  time  step 
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At  is  known  from  Equation  (10).  Equation  (42)  or  Equation 
(44)  yields  the  increase  in  the  porewater  pressure  for  the 
time  interval  At.  The  net  effective  vertical  stress  is  modi¬ 
fied  using  the  value  obtained  from  Equation  (42)  for  in¬ 
crease  of  porewater  pressure.  Equation  (15)  can  then  be  used 
to  get  the  new  value  of  shear  modulus. 


2.9  ALGORITHM  FOR  EFFECTIVE  STRESS  ANALYSIS 


Using  the  procedure  outlined  for  solving  dynamic  equa¬ 
tion  of  motion  in  conjunction  with  the  diffusion  equation, 
the  algorithm  for  the  computer  program  implementing  the  'en¬ 
gineering  approach'  is  as  follows: 

1.  Set  up  the  mass  matrix  using  Equati on ( 17 ) . 

2.  Set  up  initial  stiffness  matrix  using  Equation(l) 
through  Equation(3)  and  Equation(18)  through  Equation 
(21). 

3.  Set  up  damping  matrix  using  Equat i on  (  25  ) . 

4.  Set  up  the  load  vector  using  Equat i on ( 22 ) . 

5.  For  suitable  time  interval  use  Newmark's  method  to 
find  displacements  of  the  lumped  masses.  Also  calcu¬ 
late  shear  strains  using  these  displacements  for  each 
layer  of  the  stratum. 

6.  Repeat  step  5  till  stage  for  changing  shear  modulus 
is  reached.  Find  the  maximum  shear  strain  upto  this 
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7.  Use  Equation  (10)  find  the  change  in  the  volumetric 
strain  for  the  largest  value  of  shear  strain  and  add 
it  to  the  initial  volumetric  strain  to  get  the  cumu¬ 
lative  volumetric  strain. 

8.  Using  the  value  of  change  in  volumetric  strain  obtain 
the  distribution  of  porewater  pressure  at  that  stage. 

9.  Obtain  the  new  values  of  net  effective  porewater 
pressure  and  from  Figure  (2)  and  Equation  (15)  find 
the  new  values  of  shear  modulus  for  each  layer. 

10.  Update  stiffness  and  hence  damping  matrix  for  new  va¬ 
lues  of  shear  modulus. 

11.  Repeat  steps  1  through  10. 

The  computer  program  based  on  the  above  scheme  is  de¬ 
scribed  in  the  Appendices. 


SECTION  III 


EXAMPLES  OF  COMPUTER  PROGRAM  USAGE 

3.1  ILLUSTRATIVE  PROBLEM 

3.1.1  Introduction 

The  computer  program  'SAND'  described  in  the  appendic¬ 
es  was  used  to  solve  a  problem  previously  analysed  by  Finn 
[6].  A  saturated  layered  sand  deposit  is  shown  in  Figure 
(6).  The  50  ft  deep  sand  stratum  had  relative  density  of  45% 
with  water  table  5  ft  below  the  surface.  Following  [6],  the 
basic  properties  of  sand  used  were: 


1.  Angle  of  shear  resistance  -  40° 

2.  Unit  weight  of  saturated  sand  - 122.0  pcf 

3.  Buoyant  weight  of  sand  - 60.0  pcf 

4.  Specific  Gravity  of  sand - -  2.65 


The  following  relations  were  used  to  calculate  void  ratio 
and  dry  density  of  sand. 

e  =  ^  ^w  ^sat^^ysat  "  YJ 

Yd  ry  =  (ysat  "  ^w^  ^  ^  (**-1) 
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Figure  6:  Layered  Sand  System  for  Analysis 


where  e  =  void  ratio,  G  =  specific  gravity,  y  =  unit  wei- 
ght  of  saturated  sand,  ydr^  *  unit  weight  of  dry  sand  and  yw 
=  unit  weight  of  water.  In  addition,  following  numerical  va¬ 
lues  used  by  Finn  [6]  were  taken  for  analysis  purpose. 

1.  Total  time  of  excitation  -  10  seconds 

2.  Time  interval  (changing  shear  modulus)  0.5  seconds 

3.  Time  interval  (step-by-step  integration)  0.01  seconds 


4.  Time  interval  (diffusion  equation) 


0.5  seconds 


5.  Number  of  sand  layers 


6.  Number  of  spatial  divisions 


7.  Stiffness  proportional  damping  was  assumed. 

8.  Damping  ratio  (defined  as  ratio  of  damping  to  criti¬ 
cal  damping)  :  Several  different  values  were  used. 
These  included  0.02,  0.10  and  0.15. 


3.1.2  Exampl e  1 


A  sinusoidal  acceleration  with  a  frequency  of  2  cycles 
per  second  and  with  a  maximum  acceleration  of  0.805  ft/s  was 
applied  at  the  base  of  the  soil  system  (Figure  7).  For  this 
case,  in  addition  to  the  three  different  damping  ratios 
listed  in  item  8  of  3.1.1,  a  damping  coefficient  indepedent 
of  frequency  was  used.  This  corresponds  to  a  value  of  0.02 
for  the  constant  a  jn  Equation  (25). 


3 


3.1.3  Examp  1 e  2 

The  sand  stratum  described  above  was  subjected  to  the 
earthquake  excitation  of  the  first  10  s  of  the  North-South 
component  of  the  El  Centro  earthquake  (1940)  scaled  to  O.lg 
maximum  acceleration  (Figure  8). 

3.2  DISCRETIZED  MODEL 

The  dynamic  response  was  obtained  for  the  sand  system, 
using  lumped  mass  representation  shown  in  Figure  (9).  The 
ten  sand  layers  were  lumped  into  ten  masses  using  the  ap¬ 
proach  illustrated  in  Figure  (3).  The  masses  were  numbered 
,  beginning  with  1  at  the  top,  to  10  at  the  intersection  of 
9  th  and  10  th  layer.  The  ground  water  was  assumed  to  be  at 
the  top  of  the  second  layer  i.  e.  the  soil  system  was  satu¬ 
rated  from  5  ft.  below  the  surface.  For  spatial  division  in 
diffusion  equation,  the  saturated  part  of  the  system  was  di¬ 
vided  into  18  equal  layers  each  of  2.5  ft.  The  nodes  for  the 
finite  difference  scheme  were  numbered  beginning  with  1  at 
the  base  to  19  at  the  water  table  level. 

3.3  METHOD  OF  SOLUTION 

Algorithm  established  in  Section  II  was  used  for  ana¬ 
lysing  the  dynamic  response  of  the  soil  system.  Newmark's 


method  with  5  =  0.5  and  a=  0.25  was  used  in  solution  of 
dynamic  equation  while  explicit  and  implicit  methods  were 
applied  for  solving  the  diffusion  equation.  After  liquefac¬ 
tion  takes  place  at  any  layer  in  the  system  no  shear  stress¬ 
es  can  be  transmitted  across  that  layer.  However,  in  the 
computer  program,  to  control  the  computational  stability, 
residual  stiffness  was  retained  by  specifying  that  the  ver¬ 
tical  effective  stress  shall  not  be  less  than  10"^  of  the 
initial  effective  stress  4'  [61. 


SECTION  IV 


RESULTS  OF  ANALYSIS 


4.1  EXAMPLE-1 

Both  the  explicit  and  the  implicit  method  were  used 
for  analysis  of  dynamic  response  of  the  system  in  Example-1. 
The  results  were  obtained  for  different  values  of  damping 
ratio  viz.  0.02,  0.10  and  0.15  for  the  explicit  method  and 
0.02  and  0.10  for  the  implicit  method.  Furthur,  the  problem 
was  also  studied  using  the  damping  coefficient  to  be  equal 
to  0.02,  thus  making  it  f requency- i ndependent . 

4.1.1  Surface  Acceleration 

The  acceleration  response  for  the  entire  excitation 
period  was  obtained  using  the  explicit  method  and  the  impli¬ 
cit  method  for  the  case  £=  0.02.  Graphical  rep  res ent at i on 
of  this  for  explicit  method  is  shown  in  Figure  (10)  and  that 
for  the  implicit  method  in  Figure  (11).  The  values  of  the 
maximum  surface  acceleration  for  the  various  cases  are  shown 
in  Table  ( 1 ) . 
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4.1.2  Shear  Stress  Variation 


The  shear  stress  variation  for  the  layer  that  lique¬ 
fied  first  is  shown  in  Figures  (12)  and  (13)  for  the  expli¬ 
cit  and  the  implicit  method  respectively.  This  is  for  the 
case  of  ^  *  0.02.  The  maximum  values  of  the  shear  stress  for 
the  liquefied  layer  for  different  values  of  damping  ratios 
a  re  shown  in  Table  ( 1 ) . 


4.1.3  Pore  Pressure  Distribution 


The  effect  of  the  pore  pressure  build-up  on  the  effec¬ 
tive  stress  of  the  system  for  different  values  of  damping 
ratio  and  for  constant  damping  coefficient  is  shown  in  Fig¬ 
ures  (14)  through  (20).  Table  (2)  shows  the  tine  and  loca¬ 
tion  of  liquefaction  for  various  cases  studied.  Results  for 
pressure  distribution  are  listed  in  Tables  (3)  through  (9). 


4.2  EXAMPLE-2 


A  similar  analysis  was  carried  out  for  Example-2  but 
here  only  one  value  of  damping  ratio  equal  to  0.02  was  used. 
The  variation  for  the  surface  acceleration  for  the  explicit 
method  and  the  implicit  method  are  shown  in  Figure  (21)  and 
(22). 
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Figure  14:  Example-1  Pressure 
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Figure  15:  Example-2  Pressure  Distribution  (  £ 
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Figure  17:  Example-1  Pressure  distribution  (  £  =  0.10) 
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Figure  18:  Example-1  Pressure  Distribution  (a. 
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Figure  20:  Example-1  Pressure  Distribution  (  £=  0.15) 


TABLE  2 


Example-1  Time  and  Location  of  Liquefaction 


Damping  ratio 
or 

coef f i ci ent 

Explicit  Method 

Implicit  Method 

Layer 

T  i  me 
(sec ) 

Layer  Time 

(sec) 

OsJ 

O 

• 

o 

II 

5 ,6  ,7 ,8 

3.0 

£  =  0.10 

6 

4.0 

No  liquefaction 

=  0.02 

5,6 

5.0 

-do- 

£  =  0.15 

5,6 

6.5 

TABLE  3 


5.0 

7.5 

2.5 

7.5 


22.5 

27.5 

32.5 

37.5 

42.5 

47.5 


Example-1  Pressure  Distribution  (£ 
( Expl i ci t  Method  ) 


Pore  Pressure  at  Time 


Depth 

Initial 

Effective  Stress 

( feet ) 

(psf ) 

239.30 

627.60 

925.60 

1223.60 

1521.60 

1819.60 

2117.6 

2415.6 

2713.6 

3011.6 


636.6 

948.0 

1209.6 

799.2 

1211.5 

1633.4 

943.0 

1447.0 

2021.2 

1065.1 

1642.5 

2303.8 

1163.5 

1789.4 

2452.7 

1237.7 

1884.5 

2492.0 

1288.0 

1929.2 

2458.9 

TABLE  4 


Example-1  Pressure  Distribution  (  <f  =  0.02) 


( Impl i ci t  Method  ) 


Depth 

Initial 

Effective  Stress 

Pore  Pressure 

at  Time 

( feet ) 

(psf ) 

1 

jj|j|S|y 

3.0 
(sec  ) 

4.0 
(sec  ) 

5.0 
(sec  ) 

5.0 

239.30 

0.0 

0.0 

0.0 

0.0 

m 

627.60 

116.8 

jjjgjj 

336.3 

12.5 

925.60 

90.7 

200.8 

336.1 

529.3 

17.5 

1223.60 

124.1 

277.9 

473.5 

795.5 

1008.0 

22.5 

1521.60 

154.3 

348.3 

602.7 

1092.7 

1507.4 

27.5 

1819.60 

181.3 

411.0 

719.9 

1391.0 

1856.9 

32.5 

2117.60 

204.9 

465.3 

820.7 

1628.5 

2083.5 

37.5 

2415.60 

225.1 

510.4 

901.5 

1772.4 

2233.0 

42.5 

2713.60 

241.9 

545.9 

960.0 

1823.2 

2304.6 

47.5 

3011.60 

255.6 

566.9 

984.4 

1772.7 

2238.3 
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TABLE  5 

Example-1  Pressure  Distribution  (£=  0.10) 
(Explicit  Method ) 


Depth 

Initial 

Effective  Stress 

Pore  Pressure  at  Time 

( feet ) 

(psf ) 

1.0 
(sec ) 

1.5 
(sec  ) 

2.0 
(sec  ) 

3.0 
(sec  ) 

4.0 
(sec  ) 

5.0 

239.30 

0.0 

0.0 

0.0 

0.0 

0.0 

7.5 

627.60 

97.7 

153.8 

209.6 

341.6 

431.2 

12.5 

925.60 

167.7 

263.7 

361.2 

600.2 

771.9 

17.5 

1223.60 

229.7 

362.7 

499.7 

845.8 

1131.6 

22.5 

1521.60 

285.7 

452.4 

625.6 

1073.8 

1  1507.6 

27.5 

1819.60 

335.6 

535.1 

737.6 

1276.0 

1850.7 

‘  32.5 

2117.60 

379.7 

601.6 

834.3 

1445.0 

2093.0 

37.5 

2415.60 

417.5 

660.6 

914.8 

1575.1 

2213.2 

42.5 

2713.60 

449.2 

708.4 

978.4 

1667.3 

2253.0 

47.5 

3011.60 

474.9 

745.4 

1025.6 

1716.4 

2252.9 

Depth 


Initial 

Effective  Stres 


Pore  Pressure  at  Time 


( feet ) 

(psf ) 

I 

1 

4.0 

(sec) 

6.0 

(sec) 

MM 

Hu 

5.0 

239.30 

0.0 

0.0 

0.0 

0.0 

0.0 

627.60 

HQQ 

197.7 

308.2 

1 

452.0 

12.5 

925.60 

167.7 

346.0 

553.9 

671.1 

781.5 

■Hy| 

1223.60 

231.7 

485.3 

803.4 

965.0 

1102.4 

|8H 

1521.60 

290.0 

615.2 

1070.4 

1277.6 

1433.3 

27.5 

1819.60 

342.2 

732.7 

1355.1 

1619.4 

1783.7 

32.5 

2117.6 

387.9 

834.0 

1596.0 

1907.2 

2093.8 

3  7.5 

2415.6 

426.6 

916.1 

1727.5 

2039.8 

Sv.-.v" 

42.5 

2713.6 

458.1 

gUfg 

1748.8 

2078.1 

2408.0 

47.5 

3011.60 

477.9 

1002.3 

1684.1 

2020.6 

2411.1 

Initial 

Effective  Stress 


Pore  Pressure  at  Time 


The  shear  stress  variation  for  the  layer  that 
liquefied  first  is  shown  in  Figure  (23)  for  the  Explicit 
method  and  in  Figure  (24)  for  the  implicit  method.  The  pore 
pressure  distribution  is  shown  in  Figures  (25)  and  (26)  for 
the  explicit  and  the  implicit  method  respectively.  Tables 
(10)  and  (11)  list  the  pore  pressure  values. 


4.3  SOME  OBSERVATIONS 

4.3.1  Example-1 

The  variation  of  surface  acceleration  over  the  period 
of  excitation  shown  in  Figures  (10)  and  (11)  indicated  that 
the  surface  acceleration  gradually  built  up  and  reached  its 
maximum.  After  this,  it  dropped  rapidly  till  liquefaction, 
when  it  reached  near  zero  values.  The  shear  stress  variation 
in  layer  No.  6,  which  liquefied  first,  showed  that  the  shear 
stress  also  increased  initially  till  it  reached  its  peak  and 
then  reduced  sharply  and  attained  near  zero  values  at  lique¬ 
faction.  Figures  (14)  through  (20)  indicate  that  the  pore 
pressure  increased  slowly  in  the  initial  stages  of  excita¬ 
tion.  However ,  a  stage  is  reached  at  which  its  value  became 
equal  to  that  of  the  initial  effective  stress.  That  point 
gives  the  liquefaction  time  and  location. 


TABLE  10 

Example-2  Pressure  Distribution  (  £  =  0.02} 


Depth 


12.5 


17.5 


22.5 


27.5 


32.5 


37.5 


42.5 


47.5 


( Expl i ci t  Method  ) 


Initial 


Pore  Pressure  at 


(feet)  Effective  Stress 

(psf ) 


239.30 


627.60 


925.60 


1223.60 


1521.60 


1819.60 


2117.60 


2415.60 


2713.60 


3011.60 


1.0  2.0  3.0 

( sec  )  ( sec  )  (sec) 


43.29  276.1 


72.79  487.0 


121.40  864.7 


168.6  1304.9 


184.80  1529.0 


611.2 


1114.0 


99.00  684.4  1623.4 


2100.8 


140.40  1022.7  2493.1 


156.40  1157.0  2772.9 


3046.5 


177.50  1428.2  3278.6 


3460.6 


Depth 
( feet ) 


Initial 


Pore  Pressure  at 


5.0 

7.5 

12.5 

17.5 

22.5 

27.5 

32.5 

37.5 

42.5 


E  f f ect i ve 
(psf ) 


239.30 


627.60 


925.60 


1223.60 


1521.60 


1819.60 


2117.60 


2415.60 


2713.60 


Stress 


1.0 

(sec ) 


0.0 


21.60 


36.40 


49.50 


60.60 


70.10 


78.10 


84.20 


88.70 


2.0 

(sec ) 


0.0 


133.9 


236.1 


331.6 


419.2 


512.2 


598.8 


673.4 


735.6 


3.0 
( sec 


0.0 


378.9 


786.3 


1243.1 


1565.6 


1787.4 


2026.8 


2270.5 


2413.0 


5 


3011.60 


92.00 


781.8 


2438.2 


4.3.2  Example-2 


The  acceleration  variation  in  Figures  (21)  and  (22) 
show  that  surface  acceleration  reached  values  of  more  than 
10.0  ft/second  .  The  shear  stress  variation  in  Figures  (23) 
and  (24)  shows  that  the  maximum  values  of  shear  stress  were 
near  about  1000.0  psf.  The  pressure  variation  in  Figures 
(25)  and  (26)  indicates  the  effect  of  pore  pressure  build¬ 
up.  The  liquefaction  occured  at  2.5  and  3.0  seconds  respec¬ 
tively  in  the  explicit  and  the  implicit  method. 

4.4  Influence  of  Damping  Ratio 

As  seen  in  Figures  (14)  through  (20),  increase  in  the 
damping  ratio,  delayed  the  liquefaction.  In  fact,  the  time 
for  liquefaction  for  £  =  0.15  is  more  than  twice  that  of  the 
time  with  £=  0.02.  The  values  of  the  maximum  acceleration 
and  shear  stress  also  decreased  with  the  increase  in  the 
damping  ratio.  Thus,  choice  of  damping  ratio  influenced  the 
overall  response  of  the  system  to  dynamic  excitation  signi¬ 
ficantly. 

4.5  CHOICE  OF  DAMPING  COEFFICIENT  {a  ) 

In  Example-1,  the  stiffness  proportional  damping  was 
used  by  taking  the  value  of  the  damping  coefficient  to  be 


0.02  for  the  entire  period  of  excitation.  This  was  felt  ne¬ 
cessary  to  check  the  influence  of  the  value  of  the  natural 
frequency  over  the  results.  As  was  observed,  while  the  time 
for  liquefaction  increased  considerably  for  the  explicit 
method,  in  the  implicit  method  of  analysis  liquefaction  did 
not  occur  at  all.  The  values  of  the  maximum  acceleration  and 
shear  stress  were  found  to  lie  between  those  given  by  the 
analysis  with  £  *  0.10  and  £  =  0.15. 


4.6  COMPARISON  OF  RESULTS  WITH  FINN  [6] 


Finn  [6]  presented  results  of  his  analysis  for  the 
case  of  damping  ratio  equal  to  0.02.  For  this  case,  the  re¬ 
sults  from  the  implicit  method,  which  is  expected  to  give 
more  accurate  results  than  the  explicit  method,  were  com¬ 
pared  with  those  given  by  Finn  [6],  Table  (12)  shows  the  va¬ 
lues  of  the  maximum  acceleration,  maximum  shear  stress,  time 
of  liquefaction  and  the  depth  at  which  the  liquefaction  oc- 
cured  first  obtained  using  the  implicit  method  in  comparison 
with  Finn's  [6]  results. 


It  is  noted  that  the  time  of  liquefaction  differs 
widely  from  that  of  Finn  [6],  but  the  depth  of  liquefaction 
is  the  same  in  both  the  examples.  Also,  the  values  of  the 
maximum  acceleration  and  shear  stress  are  significantly 
d i f f eren  t . 


TABLE  12 


Comparison  of  Results  with  Finn  [6] 


Implicit  Method 

Finn's  Analysi  s 

Example  1 

Example  2 

Example  1 

Example  2 

Maxi  mum 

Acce 1 erat i on 

5.0 

5.2 

3.0 

13.0 

Maximum 

Shear  Stress 

500.0 

1020.0 

250.0 

400.0 

Time  of 
Liquefaction 

5.0 

3.0 

6.0 

5.0 

Depth 

27.5 

22.5 

27.5 

22.5 

*  Units  used: 


Acceleration 
Shear  Stress 
T  i  me 
Depth 


ft/second2 

psf 

seconds 
ft . 


i 


I 


I 


I 


SECTION  V 


CONCLUSIONS 

Numerical  examples  using  computer  program  "SAND"  based 
on  the  'engineering  approach'  indicate  the  following: 

1.  The  explicit  method  and  implicit  method  give  diffe¬ 
rent  results  for  the  same  time  step.  However,  use  of 
smaller  time  step  for  explicit  to  establish  conver¬ 
gence  to  results  obtained  using  the  implicit  method 
could  not  be  done,  since  the  choice  of  time  step  is 
dependent  upon  the  frequency  of  input  excitation  and 
cannot  be  arbitrarily  reduced. 

2.  The  comparison  of  results  from  the  implicit  method 
with  those  reported  by  Finn  [6]  shows  a  wide  differ¬ 
ence  in  the  values  of  maximum  surface  acceleration 
and  shear  stress.  This  could  be  due  to  one  or  more  of 
the  following  reasons: 

a)  Finn  [6]  has  not  indicated  how  the  natural  fre- 
que  cy  of  the  system  was  established.  This  is  es¬ 
sential  to  calculate  the  constant  a  in  Equation 
(24).  Since  for  stiffness  proportional  damping, 
the  damping  ratio  is  directly  proportional  to  the 
natural  frequency,  this  factor  could  be  quite  sig¬ 
nificant. 
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b)  The  expression  for  updated  values  of  the  shear  mo¬ 
dulus  was  not  same  as  that  used  by  Finn  [6].  In 
this  study,  shear  modulus  is  calculated  using 
Equation  (15)  while  Finn  [6]  had  based  it  on  Equa¬ 
tions  (4),  (6)  and  (7).  The  values  of  constants 
used  in  these  equations  are  not  available. 

c)  The  method  of  analysis  for  solution  of  diffusion 
equation  is  not  mentioned  by  Finn  [6].  In  this 
study,  both  the  explicit  and  the  implicit  methods 
were  used.  This  part  of  establishing  values  of 
pore  pressures  is  very  important  in  this  type  of 
analysis,  since  the  updated  values  of  effective 
stresses  are  based  on  pore  pressure  distribution. 

This  study  brought  forth  following  drawbacks  of  the 
engineering  approach'  viz. 

1.  The  method  is  based  on  simplifying  assumptions  such 
as  the  sequence  of  events  listed  in  SECTION  II,  which 
is  only  an  approximate  representation  of  the  physical 
phenomenon.  In  reality,  the  variation  of  pore  pres¬ 
sure  and  soil  properties  is  more  likely  to  occur  sim¬ 
ultaneously  and  a  correct  theory  must  simulate  that. 
For  an  approximate  theory  to  be  acceptable,  its  lim¬ 
its  and  validity  must  be  carefully  determined. 

2.  The  method  requires  a  large  amount  of  input  informa¬ 
tion  in  terms  of  soil  properties  and  has  several  re- 


lations  based  on  experimental  analysis.  Without  such 
information,  the  method  cannot  possibly  predict  the 
time  and  depth  of  liquefacton  of  soil  stratum.  These 
data  are  site  specific  and  would  be  expensive  to  get. 
The  pore  pressure  is  established  using  Equation  (10) 
for  calculating  change  in  volumetric  strain  and  Equa¬ 
tion  (14).  Equation  (10)  is  however  based  on  experi¬ 
mental  data.  Since  this  is  an  important  factor  to¬ 
wards  establishing  soil  properties  in  successive 
stages  of  excitation,  this  requires  to  be  carefully 
studied. 

The  method  is  established  for  the  case  of  horizontal¬ 
ly  layered  soil  strata  only  and  the  properties  are 
assumed  to  vary  only  in  the  vertical  direction.  This 
places  a  serious  restriction  on  the  use  of  this  meth¬ 
od  . 

The  soil  properties  are  modified  by  updating  the  va¬ 
lues  of  shear  modulus.  The  basis  for  this  updating 
needs  to  be  examined,  since  this  is  related  to  the 
cost  of  computer  program.  The  method  is  found  to  be 
expensive  from  this  point  of  view. 

The  study  established  that  the  results  of  the  program 
are  dependent  on  the  method  of  analysis.  Also  the 
time  step  for  solution  of  diffusion  equation  is  based 
on  the  frequency  of  the  input  excitation.  Since  this 
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is  related  to  establishing  other  soil  properties,  it 
is  not  possible  to  set  up  a  sequence  of  solutions 
with  systematically  decreasing  time  step  size,  to  es¬ 
tablish  convergence  to  a  unique  result. 


section: vi 
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Appendix  A 


PROGRAM  STRUCTURE  AND  ORGANIZATION 

A . 1  INTRODUCTION 


The  computer  program  consists  of  a  main  unit  and  sev¬ 
eral  subroutines.  Figure  (A.l)  shows  the  nesting  of  subrou- 
ti nes . 

A . 2  MAIN 

In  this  unit,  the  control  information  is  read  in.  The 
dynamic  storage  allocation  is  set  up  and  locations  of  vari¬ 
ous  arrays  specified. 

A. 3  I NPT 

This  subroutine  is  called  by  the  unit  MAIN.  In  this 
subroutine,  the  dynamic  storage  allocation  is  completed  and 
storage  requirement  for  problem  under  consideration  is 
checked.  This  subroutine  also  forms  the  main  body  of  the 
program  and  calls  in  various  subroutines  for  specific  opera¬ 
tions. 
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DYNR 


REDUCE 


J 


RMRX 


SOLVE 


A. 4  SPLINE 


This  subroutine  is  called  by  the  unit  I  MPT.  In  this 
subroutine,  the  spline  interpolation  of  the  data  points  is 
carried  out. 

A. 5  TLOAD 


This  subroutine,  called  by  the  unit  INPT,  sets  up  the 
load  vector  in  the  dynamic  equation  of  motion,  based  on  the 
input  base  excitation.  While  the  values  of  the  acceleration 
for  harmonic  excittation  are  evaluated,  the  values  for  other 
type  of  excitation  are  to  be  read  in.  The  choice  is  to  be 
indicated  by  using  ICOUNT  =  1  for  harmonic  base  input  and 
ICOUNT  =  1  for  other  types.  The  input  frequency  for  harmon¬ 
ic  excitation  is  to  be  expressed  as  a  multiple  of  = 
3.1415927  and  this  multiple  is  indicated  in  form  of  the 
vari abl e  AMO  MEG .  • 

A. 6  STIFF 


This  subroutine  is  called  by  the  unit  INPT  and  sets  up 
the  system  stiffness  matrix  for  every  cycle  of  periodic  up¬ 
dating  of  soil  properties. 


A. 7  DAMP 


This  subroutine  as  also  called  by  the  I N P T  unit  and 
sets  up  the  system  damping  matrix  for  periodic  updates. 

A. 8  DYNA 

This  subroutine  is  called  by  the  unit  INPT  for  solving 
the  dynamic  equation  of  motion  and  gives  displacements,  ve¬ 
locities  and  accelerations.  The  shear  strains  are  also  eval¬ 
uated  by  this  subroutine. 

A. 9  SOURCE 

This  subroutine  is  called  by  the  unit  INPT  and  calcu¬ 
lates  the  source  term  in  the  diffusion  equation  expressed  as 
function  of  change  of  volumetric  strain. 

A. 10  DIFF 

This  subroutine,  called  by  INPT,  solves  the  diffusion 
equation  using  both  the  explicit  and  implicit  method.  This 
is  achieved  by  using  IFLAG  =  1  for  the  explicit  method  and 
IFLAG  =  1  for  the  implicit  method. 


A. 11  SPEVAL 


This  subroutine,  called  by  the  unit  INPT,  calculates 
the  interpolated  value  using  spline  interpolation  already 
established  in  SPLINE. 

A. 12  STRESS 

This  subroutine  is  called  by  INPT  and  evaluates  the 
pore  pressure  and  effective  stresses  at  the  end  of  each  cy¬ 
cle  of  periodic  update. 

A. 13  TRI 

This  subroutine,  called  by  the  subroutine  DYNA,  per¬ 
forms  standard  Gaussian  elimination  on  the  coefficient  ma¬ 
trix  of  a  system  of  linear  equations.  No  pivoting  is  used. 
It  reduces  the  matrix  to  triangular  form. 

A. 14  REDUCE 

This  subroutine  is  called  by  the  subroutine  DYNA  for 
row  operations  on  the  vector  in  a  system  of  linear  equations 
and  carries  out  back-substitution  to  evaluate  the  unknown 
di spl acements . 


AD-A150  769  ft  COMPUTER  PROGRAM  FOR  DYNAMIC  RESPONSE  OF  LAYERED  2/2 

SATURATED  SANDS.  (U)  OHIO  STATE  UNIV  RESEARCH 
FOUNDATION  COLUMBUS  M  S  HIREMATH  ET  AL.  JUN  84 
UNCLASSIFIED  OSURF-715107-84-3  AFOSR-TR-85-0088  F/G  9/2  NL 


This  subroutine  is  called  by  the  subroutine  QYNA  to 
calculate  the  maximum  shear  strain  for  the  cycle  of  period¬ 


ic  update  of  soil  properties. 

A. 16  SOLVE 

This  subroutine,  called  by  the  subroutine  DIFF,  solves 
tri-diagonal  system  of  equations. 


Appendix  B 
PROGRAM  USAGE 


B.l  PROGRAM  IDENTIFICATION 


Name : 

Programmer : 


SAND 


Date : 


M.  S.  Hiremath,  Ohio  State  University, 
Columbus,  Ohio. 

April  1984. 


B.2  INPUT  DATA 


Input  to  the  program  is  a  sequence  of  punched  cards  in 
the  following  order  and  format. 


B.2.1  Problem  Identification  (18A4) 


One  card  carries  a  descriptive  title  to  ident 
ify  the  job. 


B.2. 2  Control  Information 


These  sets  of  cards  will  provide  the  following 
control  information  for  the  problem  being  solved. 
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i.  First  Card  (4110) 

Informati on 

Total  number  of  layers  (NO)  . 

Number  of  dry  layers  (LAYDRY)  . 

Number  of  data  points  (NN)  . 

(Used  in  i nterpol at i on  of  curve  in 
Figure  ( 2 ) ) . 

Number  of  spatial  divisions  (N)  . 

ii.  Second  Card  (8F10.0) 

Informati  on 

Depth  of  sand  system  (Z)  . 

Acceleration  due  to  gravity  (EG)  .... 

Damping  rati o  (ZETA)  . 

Unit  weight  of  saturated  sand  (GAMSAT) 

Unit  weight  of  water  (GAMW)  . 

Specific  gravity  of  sand  (SPGRA)  .... 

Angle. of  shear  resistance  (PHI)  . 

Relative  density  (RELDEN)  . 

iii.  Third  Card  (10F5.0) 

Informati on 

Values  of  shear  strain  in  Figure  (2) 
(In  all  NN  values  are  to  be  provided 
one  value  for  in  every  five  columns) 
Array  is  denoted  by  (X). 


Col  umns 
1  -  10 
11  -  20 
21  -  30 

31  -  40 


Col umns 
1  -  10 
11  -  20 
21  -  30 
31  -  40 
41  -  50 
51  -  60 
61  -  70 
71  -  80 


Col umns 
1  -  50 
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iv.  Fourth  Card  (10f5.0) 

I nf ormat i on  Col umn  s 

Values  of  K  in  Figure  (2)  corresponding  1-50 
to  the  values  in  array  (X).  NN  values 
are  to  be  provided,  one  in  every  five 
columns.  The  Array  is  denoted  by  (V) 

v.  Fifth  Card  (3F10.0) 


I nf ormat i on  Col umns 

Value  of  'm‘  in  Equation  (11)  (AM)...  1  -  10 

Value  of  'n'  in  Equation  (11)  (AN)...  11  -  20 

Value  of  'k'  in  Equation  (11)  (AK)...  21  -  30 

vi .  Sixth  Card  (F10.0) 

Information  Col umns 

Initial  shear  strain  (SSTI)  .  1-10 


vii.  Seventh  Card  (4F10.0) 


I nf ormat i on  Col umns 

Permeability  of  sand  system  (PERM)  ..  1-10 

Depth  of  saturated  portion  (XL)  .  11  -  20 

Time  step  for  updating  properties  (D)  21  -  30 


Total  time  period  of  excitation  (TMAX)  31  -  40 


vi i i .  Ei ghth  Card  ( i 10  ) 

In  forma t i on  Col umns 

Number  of  divisions  of  TMAX  .  1-10 


i x .  Ninth  Card  ( 3F 10 .0  ) 

I n  format i on 

Multiple  of  1  n  '  in  input  frequency 
(ANOMEG) 

Time  step  for  calculating  load  (DELT). 
Maximum  Acceleration  (AA)  . 

x.  Tenth  Card  (4F10.0) 


Information  Col umns 

Time  step  for  Newmark's  method  (TSTEP)  1-10 

Newmark's  constant  ($)  (DELTA)  .  11  -  20 

Newmark's  constant  (a  )  (ALPHA)  .  21  -  30 


Time  period  for  Newmark’s  method  (TFIN)31  -  40 

xi .  El eventh  Card  (110) 


I nf ormat i on  Col umns 

Number  of  divisions  of  TFIN  .  1-10 


B.2.3  Initial  Value  Data 

Corresponding  to  each  lumped  mass,  one  card  is 


Col umns 
1  -  10 

11  -  20 
21  -  30 


to  be  provided,  giving  the  value  of  initial  displacem¬ 
ent  and  initial  velocity. 


I n  f ormat  i on 


Col umns 


Initial  Displacement  (UDI(ND))  .  1  -  10 

Initial  velocity  (UVI(ND))  .  11  -  20 

B  .  2.4  Soi  1  Constants 

Soil  constants  in  Equation  (10)  are  to  be  sup¬ 
plied  as  foil ows . 

I nf ormat i on  Col umns 

Value  of  (Cl)  .  1  -  10 

Value  of  C2  (C2)  .  11  -  20 

Value  of  C3  (C3)  . 21  -  30 

Value  of  C4  (C4)  .  31  -  40 

B.2.5  Choice  of  the  Method 

Choice  of  the  method  and  type  of  base  input  is 
to  be  indicated  using  following  two  cards. 

i .  First  Card  (110) 

I nformati on  Col umns 

IFLAG  =  1  for  Explicit  method  1-10 

»  2  for  Implicit  method 


i i .  Second  Ca rd  (110) 


Appendix  C 

DICTIONARY  OF  VARIABLE  NAMES 


Following  is  the  list  of  variables  used  in  the  compu¬ 
ter  program  SAND.  The  dimensions  of  the  arrays  are  indicated 
and  may  be  altered  to  suite  the  problem  to  be  analysed. 


VARIABLE 


DESCRIPTION 


AKO 

Coefficient  of  earth  pressure  at 

rest . 

AK2 

Interpolated  value  of  the  coeffcient 

Y 

A,  B,  C 

N 

Diagonals  in  a  tri-diagonal  system. 

A1 

Stiffness  damping  coeffcient. 

A ( ND , ND  ) 

Coefficient  matrix  in  a  set  of 

1 i near  equati ons . 

Aj ,  thr  A i 

Constants  in  NEWMARK'S  method. 

B  ( ND  ) 

The  known  vector  associated  with 

a  set  of  linear  equations. 

BETA(N) 

Constant  in  explicit  method. 

BET(ND) 

Value  of  the  constant  term  in  the 

explicit  method  associated  with  each 

sand  layer. 
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CHVD(ND) 


Change  in  the  volumetric  strain 
during  the  cycle  under  consideration. 


DEN(NO)  Denominator  in  the  expression  for  E 

DELVD(ND)  Change  in  volumetric  st ra i n ( =CH VD  ) 

EE  Void  ratio 

ER(ND)  One-dimensional  rebound  modulus. 

EER ( N )  Value  of  E  associated  with  each 

source  term  in  the  diffusion  equation 
EPSVD(ND)  Cumulative  volumetric  strain  for 

each  layer. 

ESTI(ND)  Effective  stress  for  each  layer. 

EKI(ND)  Value  of  the  spring  constant  for  each 

1 ayer . 

Interpolated  value  of  K^. 

Array  for  load  vector  in  dynamic 
equation  of  motion  for  entire  period 
of  excitation. 

Values  of  the  pore  pressure  in  the 
diffusion  equation. 

Array  containing  second  derivatives  i 
interpolating  cubic  spline. 

GAMDRY  Unit  weight  of  dry  sand. 

GAMA ( NO )  Unit  weight  of  each  sand  layer. 

GMO(ND)  Initial  value  of  the  shear  modulus 


F 

FF ( ND ,NT+1 ) 

FR(ND) 

FDP(NN) 


61  (NO) 


G  ( NO ) 

H  ( ND ) 

KKK 

PPI  (ND) 

Q(N) 

R 

SM ( ND , ND  ) 

SK(ND.ND) 

SC(ND.ND) 

SIGZER(ND) 

SST(ND) 

SIGMAO(ND) 

SIGMAV(ND) 

SHE AR ( NO  ) 
SKK(ND) 


for  each  1 ayer . 

Values  of  the  shear  modulus  for 
each  layer. 

Duplicate  of  GI(ND). 

Depth  of  each  layer. 

Counter  for  number  of  cycles. 

Pore  pressure  at  the  center  of  each 
1 ayer . 

Source  term  in  the  diffusion  equation. 

Right  hand  side  vector  in  tri-diagona! 
system  in  Subroutine  S P L I N E > 

Mass  matrix  in  Dynamic  analysis. 
Stiffness  matrix  in  Dynamic  analysis. 
Damping  matrix  in  Dynamic  analysis. 
Initial  values  of  the  shear  stresses. 
Values  of  the  shear  strains. 

Mean  normal  effective  stress  for  each 
1 ayer . 

Net  effective  stress  at  the  end  of  each 
cycl e . 

Shear  stress  in  each  layer. 

Dupl i cate  for  E  K I ( ND  ) 


TSTI (ND) 


Initial  total  stress  for  each  layer. 


TMO ( NO  ) 


Value  of  the  initial  maximum  shear 


T(N-l) 

UAI(ND) 

UD(ND) 

UV(ND) 

UA(ND) 

WN 

XX 

XL 

Y  ( NO ) 

Z 

ZETA 


stress  for  each  layer. 

Temporary  array  in  subroutine  DIFF 

Initial  accel erati on  . 

Displacements  of  masses. 

Velocities  of  masses. 

Accelerations  of  masses. 

Natural  frequency  of  the  stratum. 

Value  of  the  maximum  shear  strain  for 
which  interpolated  value  of  K g  is 
calculated  in  subroutine  SPEVAL. 
Length  of  spatial  domain  in 
exp? icit  method. 

Data  points  for  cubic  spline 
interpolation  used  for  finding  . 

Total  depth  of  sand  stratum. 

Damping  ratio. 


o  o  o  o  o  o  o 


Appendix  D 
PROGRAM  LISTING 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

COMPUTER  PROGRAM  "SAND" 


FOR  DYNAMIC  RESPONSE  OF  FLUID  SATURATED  LAYERED 
SAND  STRATUM 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

IMPLICIT  REAL*8(A-H  ,0-Z  ) 

COMMON  RA ( 18000  ) 

COMMON  /AAA/  ND  , L A YD R Y  , NN  , N  , NM 
COMMON  / B B B /  NT  ,  ND  I V  ,  NTP  ,  ND  I  VP 
COMMON  /CCC/  N36 
C 

DIMENSION  HEAD (18 ) 

C 

READ  (5,20)  HEAD 
C 

READ  (5,30)  ND,LAYDRY,NN,NLD 
C 

READ  (5,40)  NT  ,  ND  I V 
C 

N  =  NLD*( ND-LAYDRY  ) 

NM  =  N- 1 

NTP  =  NT  +  1 

NDIVP=NDIV+1 

C 

WRITE  (6,50)  HEAD,ND,LAYDRY,NN,NLD,N,NT,NDIV 
C 

N 1  =  1 

N2=N I +ND 

N3=N2+ND 

N4=N3+ND 

N5=N4+ND 

N6=N5+ND 

N7=N6+ND 

N8=N7+ND 

N9=N8+ND 

N10=N9+ND 

N11=N10+ND 


99 


o  o  o  o  o  o  o 


:iTr 


N 12  =N1 1 +ND 

N13=N12+ND 

N14=N13+ND 

N15=N14+ND 

N16=N15+ND 

N 1 7 =N16+ND 

N18=N17+ND 

N19=N18+ND 

N20=N19+ND 

N21 =N20+NN 

N22=N21+NN 

N23=N22+NN 

N24=N23+ND 

N25=N24+ND*ND 

N26=N25+ND*ND 

N27=N26+ND*ND 

N28=N27+ND 

N29=N28+ND 

N30=N29+ND 

N31 =N30+ND*NTP 

N32=N31+N 

N33=N32+N 

N34=N33+N 

N35=N34+N 

N36=N35+ND 

CALL  INPT  (RA(N1),RA(N2),RA(N3),RA(N4),RA(N5),RA(N6),RA(N7),RA(N8) 
1,RA(N9),RA(N10),RA(N11),RA(N12),RA(N13),RA(N14),RA(N15),RA(N16),RA 
2(N17),RA(N18),RA(N19),RA(N20),RA{N21),RA(N22),RA(N23),RA(N24),RA(N 
325) ,RA(N26 ) ,RA{N27 ) ,RA(N28) ,RA(N29) ,RA(N30 ) ,RA(N31 ) ,RA(N32  )  ,RA(N33 
4),RA(N34),RA(N35)) 

STOP 


FORMAT  (18A4) 

FORMAT  (415) 

FORMAT  (215) 

FORMAT  ( 1H1 , / /18A4//5X , 1  TOTAL  NUMBER  OF  LAYERS 

1=  ' ,I5/5X, 'NUMBER  OF  DRY  LAYERS  =  * , I 5/5X , ' NU 

2MBER  OF  DATA  POINTS  FOR  INTERPOLATION  =  \I 5/5X NUMBER  OF  DIVIS 
3IONS  OF  EACH  LAYER  =  ' , 1 5/5X  , ' TOTAL  NUMBER  OF  SPATIAL  DIVI 

4SI0NS  =  '.I5/5X, 'NUMBER  OF  DIVISIONS  OF  EXCITATION  PERIOD 

5  =  ' ,I5/5X, 'NUMBER  OF  DIVISIONS  OF  EACH  CYCLE  PERIOD  =  ’,15///) 


o  n 
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SUBROUTINE  INPT  ( H  ,GAMA , TST I ,PP  I  , EST I  , DEN , E R , GMO , TMO , G I , G ,B IG , S IGM 
1AV,SIGMA0,SSTR,SHEARM,CHVD,SIGZER,EKI,X,Y,FDP ,AK2  ,SM, SK ,SC .EPSVD ,U 
2D  I  ,  UV I  ,FF,Q,FR, BETA, EER, BET) 


IMPLICIT  RE AL *8 ( A-H  , 0-Z  ) 

COMMON  RA { 18000  ) 

COMMON  /AAA/  ND .LAYDRY  ,NN  ,N  ,NM 
COMMON  /BBB /  NT , ND I V  ,  NTP  ,  ND I  VP 
COMMON  /CCC/  N36 
COMMON  /ONE/  KKK 
COMMON  /TWO/  ILIQ 

COMMON  /THREE/  DELT ,TMAX ,AA .ANOMEG 
COMMON  /FOUR/  TSTEP , DELTA , AL PHA  , TF I N 
COMMON  /FIVE/  C1,C2,C3,C4 
C 

DIMENSION  H ( ND  )  ,  GAMA(ND),  TSTI(ND),  PPI(ND),  ESTI(ND),  DEN(ND),  E 
1R(ND),  GMO(ND),  TMO(ND),  GI(ND),  G ( ND  )  ,  BIG (ND  )  ,  SIGMAV ( ND  )  ,  SIGMA 
20 ( ND ) ,  SSTR(ND),  SHEARM(ND),  CHVD(ND),  SIGZER ( ND ) ,  EKI(ND),  X(NN), 
3  Y ( NN  )  ,  FDP(NN) 

DIMENSION  AK2 ( ND  ) ,  Q  (  N  )  ,  FR ( N  )  ,  BETA(N),  EER(N),  BET(ND),  SM(ND  ,ND 
1),  SK ( ND  ,  ND  ) ,  SC ( ND  ,  ND  ) ,  EPSVD(ND),  UDI(ND),  UVI(ND),  FF(ND.NTP) 
N37=N36+NN 
N38=N37+NN 
N39=N38+NN 
N40=N39+NN 
C 

N41 =N40+NT 
C 

N42=N41+ND 

N43=N42+ND 

N44=N43+ND 

N45=N44+ND 

N46=N45+ND 

N47=N46+ND*NDIVP 

N48=N47+ND 

N49=N48+ND*ND 

N50=N49+ND 

N51=N50+ND 

N52=N51+ND 

C 

N53-N52+N 
N54=N53+N 
N55=N54+N 
N56=N55+N 
H57  =  N56  +  N 
N58=N57+N 
C 

N59=N58+ND 

C 
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WRITE  (6,300)  N59 
NT0T= 18000 

I F ( N59  .  LE . NTOT  )  GO  TO  10 
WRITE  (6,310) 

CALL  EXIT 

CONTINUE 

READ  (6,320)  Z , EG , ZETA  .GAMSAT ,GAMW , SPGRA  ,  PH  I  ,  RE LDEN 
WRITE  (6,330)  Z , EG  ,ZET A .GAMSAT , GAMW , SPGRA  ,  PHI  ,  RE LDEN 

READ  (5,340)  (X ( I ) ,  1  =  1  ,ND  ) 

WRITE  (6,350)  (X  ( I  )  ,  1  =  1  ,ND  ) 

READ  (5,360)  ( Y ( I  )  ,  I =1  ,ND  ) 

WRITE  (6,370) 

READ  (5,380)  AM , AN  ,AK 
WRITE  (6,390)  AM  ,  AN , AK 

READ  (5,400)  SSTI 
WRITE  (6,410)  SSTI 

READ  (5,420)  PERM, XL ,DT ,TMAX 
WRITE  (6,430)  P ERM , XL , DT , TMAX 

READ  (5,440)  ANOMEG ,0ELT , AA 
WRITE  (6,450)  ANOMEG ,DELT ,AA 

READ  (5,460)  TSTEP  ,  DELTA  .ALPHA  ,TF I N 
WRITE  (6,470)  TSTEP .DELTA, ALPHA, TFIN 

DO  20  1  =  1, ND 

READ  (5,480)  UD I  ( I  )  ,  UV  I  ( I  ) 

WRITE  (6,490  )  UD I  ( I  ) ,  UV  I  ( I  ) 

CONTINUE 

READ  (5,500)  C1,C2,C3,C4 
WRITE  (6,510)  C1,C2,C3,C4 

READ  (5,520)  I  FLAG 
WRITE  (6,520)  IFLAG 

READ  (5.530)  ICOUNT 
WRITE  (6,530)  ICOUNT 


CALL  SPLINE  ( NN , X , Y , FDP , RA ( N36 ) , RA ( N37 ) , RA ( N38 ) , RA ( N39  ) ) 


Ooo  OOOO  OOOO  O  O  O  o  oooo 
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C  CALCULATE  SOIL  PROPERTIES  AND  PRINT  THEM 
C 

P I =4 . *DATAN ( 1 .  DO  ) 

RAD=PHI*PI/180. 

AK0  =  1.0-DSIN( RAD  ) 

EE  =  ( SPGRA*GAMW-GAMSAT ) / ( GAMSAT -GAMW  ) 

GAMDRY= ( GAM SAT -GAMW  )*SPGRA/( SPGRA-1 .0) 

NOTE  :  THE  SAND  LAYERS  ARE  NUMBERED  FROM  TOP 


AND  =  ND 
HH=Z/AND 

DO  30  1=1, ND 
H { I ) =HH 
30  CONTINUE 

DO  40  1=1 , LAYDRY 
40  GAMA ( I )=GAMDRY 

LAYDRl =L AYDRY+1 
DO  50  I =LAYDR1  ,ND 
50  GAMA ( I )=GAMSAT 

LAYWET=ND-LAYDRY 
LAYER= ( LAYWET+1 ) /2 
LAYMID=LAYDRY+LAYER 

GAMAC=GAMA ( LAYMID ) 


WRITE  (6,540) 

WRITE  (6,550  )  Z  ,  EG  ,  ZETA  .GAMSAT  ,  GAMW , SPGRA , PHI , AKO ,EE .GAMDRY ,HH ,REL 
1DEN 

WRITE  (6,560) 

SET  UP  THE  MASS  MATRIX 

DO  60  1  =  1, ND 
DO  60  J  =  1  ,ND 
SM ( I ,J  )=0.0 
60  CONTINUE 

CALCULATE  THE  NONZERO  ELEMENTS  OF  THE  MATRIX 

SM(1,1)=GAMA(1)*H(1)*0.5/EG 

C 

DO  70  1=2, ND 

SM ( I , I )  =  { GAMA (I-1)*H(I-1) +GAMA ( I  )  *H ( I  ) ) *0 . 5 /EG 


OOO  OOOOC1  ooo  o  o  o  oooo 


70  CONTINUE 

CALCULATE  THE  TOTAL  PRESSURE,  PORE  PRESSURE  AND  EFFECTIVE  PRESSURE 
AT  THE  CENTER  OF  EACH  LAYER 

TST I  ( 1 ) =0 . 5*H ( 1 ) *GAMA ( 1 ) 

DO  80  1=2, ND 

TSTI (  I  )=TSTI (I-1)+0.5*(H(I )*GAMA ( I )  +  H( I -1 ) *GAMA ( I - 1 ) ) 

80  CONTINUE 

PPI (1)«0.0 

PPI (2)=0.5*H{2)*GAMW 
DO  90  1=3, ND 

PPI ( I )=PPI ( 1-1 )+0.5*(H( 1-1 )+H( I ))*GAMW 
90  CONTINUE 

DO  100  1  =  1, ND 

ESTI ( I ) =TST I ( I  ) -PP  I  (  I  ) 

100  CONTINUE 


DUPLICATE  THE  ARRAY  ESTI 
DO  110  1  =  1, ND 
SIGZER ( I  )  =  ESTI  ( I  ) 

110  CONTINUE 

WRITE  (6,570) 

WRITE  (6,580)  ( TSTI  ( I ) , ESTI ( I ) , PP I ( I  ) ,  I  =  1 ,  ND  ) 


CALCULATE  DENOMINATOR  IN  EXPRESSION  FOR  ER 

ANM=AN-AM 

AMI =1. -AM 

DO  120  1  =  1, ND 

DEN ( I )=AM*AK*( ESTI ( I )**ANM) 

120  CONTINUE 


CALCULATE  GMO  AND  TMO  AT  THE  CENTRE  OF  EACH  LAYER 
FAC1  =  ( ( 2.973-EE )**2 )/( 1 .+EE  ) 

FAC2  =  DSQRT( ( 1 . +2 . *AKO ) /3 . 0  ) 

C 

DO  130  1  =  1, ND 

GMO( I )  =  14760. 0*FAC1*FAC2*DSQRT{ ESTI  (I )) 

130  CONTINUE 
C 

FAC3=((1.+AK0)*0.5*DSIN(RAD)  )**2 
FAC4=( ( l.-AK0)*0.5)**2 
FAC5=FAC3-FAC4 


OOO  OO  O  OO  O  O  OOOOOOO  Ooo 
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DO  140  1  =  1, NO 

TM0( I ) =  0 SQ RT (FAC5  ) *E  STI  (  I  ) 

140  CONTINUE 
C 

CALL  TL0A0  ( NTP , FF , NT , NO , SM , I  COUNT  ,  RA ( N40  ) ) 

C  =============================================== 

C 

DO  150  1  =  1, NO 

G I ( I ) =GM0 ( I ) / ( 1 . 0+GM0 ( I )*SSTI /TM0( I  ) ) 

150  CONTINUE 

GMID  =  GI ( LAYMID  ) 

C 

00  160  1=2, N 
FR ( I  )  =0 . 0 
160  CONTINUE 

INITIALISE  THE  CUMULATIVE  VOLUMETRIC  STRAIN. 

DO  170  1  =  1, NO 
EPSVD ( I  )  =0 . 0 
170  CONTINUE 

START  OF  THE  BIG  LOOP  WHICH  INVOLVES  SOLVING  OF  DYNAMIC  EQUATION 
OF  MOTION  AND  DIFFUSION  EQUATION  AND  SUBSEQUENT  UPDATING  OF  THE 
STIFFNESS  MATRIX  ETC. 

INITIALISE  THE  COUNTERS 


T  =  0 . 0 
IL  IQ  =  0 
KKK  =  1 

180  CONTINUE 

WRITE  (6,590)  KKK 

WRITE  (6,600)  ZETA,  ILIQ 

WRITE  (6,610) 

DO  190  1  =  1, ND 
EKI ( I )=GI ( I  )/H(  I  ) 

190  CONTINUE 

CALL  STIFF  ( EK I  ,ND  , SK  ) 


DO  200  1=1  , ND 
G ( I ) =G I (  I ) 

200  CONTINUE 

CALCULATE  THE  FUNDAMENTAL  FREQUENCY  OF  THE  SIXTH  LAYER  OF  THE 
STRATA  SINCE  IT  LIES  AT  THE  MIDDLE  OF  THE  SAND  STRATUM 


106 


ZH  = ( Z-HH  )*0 . 5 
ZH=ZH**0 . 3333 
AW  =  G ( LAYMID ) / Z  H 

WN=1.7510*DSQRT ( AW*EG /GAMAC )/(1.2*(Z-HH)**0. 83333) 
CALL  DAMP  (WN,ZETA,ND,SK,SC) 


CALL  DYNA(BIG,SM,SK,SC,FF,GI  ,H,UDI  ,UVI  , 


RA ( N41 )  ,RA(N42) , RA ( N43  ) ,  RA ( N44  ) ,RA(N45) ,RA(N46) 
RA(N47),RA(N48),RA(N49),RA(N50)) 


DO  210  1  =  1, ND 

SHEARM ( I  )  =GI ( I ) *8 1 G ( I ) 

210  CONTINUE 

CALL  SOURCE  (  ND  ,KKK  , CHVD , EPSVD , B IG , RA ( N51 ) ) 


CONVERT  SHEAR  STRAINS  INTO  PERCENTAGES 

DO  220  1=1, ND 
SSTR ( I )  =  100 . *B IG ( I  ) 

IF(SSTRU).GT.l.O)  SSTR  ( I  )  =  1 .0 
IF  (  SSTR ( I ).LT. 0.0001  )  SSTR ( I ) =0 . 0001 
220  CONTINUE 


EVALUATE  THE  VALUES  OF  ELASTIC  REBOUND  MODULUS 
FOR  EACH  LAYER 

DO  230  1  =  1  , ND 

ER  ( I  )  =  ( ESTI ( I )**AM1 )/DEN ( I ) 

230  CONTINUE 

EVALUATE  THE  CONSTANT  TERM  "BETA”  FOR  EACH  LAYER 
AN  =  N 

DX  =  XL /AN 
DO  240  1  =  1  , ND 

BET ( I )  =  ER ( I )*PERM*DT/( GAMW*DX*DX ) 

240  CONTINUE 


NM1=N-1 

NND=ND 

DO  250  1=2  , N, 2 
Q ( I  )  =CHVD { NND ) 
EER ( I ) =ER ( NND ) 


BETA ( I  )  =  BET ( NND ) 

NND=NND-1 
250  CONTINUE 

COMPLETE  THE  OTHER  ELEMENTS  OF  Q,  ER,  BETA 

DO  260  I  =  3  , NM1 , 2 
Q(I)=(Q(I-1)+Q(I+1))*0.5 
EER(  I  )  =  (EER{ 1-1 )+EER( 1  +  1 ) )*0.5 
BETA(I)=(BETA(I-1)+BETA(I+1))*0.5 
260  CONTINUE 

DO  270  1=2, N 

270  Q( I  )  =  EER( I ) *Q ( I )  *100 . 0 

CALL  DIFF  (IFLAG,N,FR,Q,BETA,NM,RA(N52),RA(N53),RA(N54),RA(N55),RA 
1 ( N56 ) , RA ( N57  ) ) 


CALCULATE  THE  VALUE  OF  K2  FOR  EACH  LAYER  BY  INTERPOLATION 


DO  280  1=1, ND 
XX  =  SSTR (  I  ) 

CALL  SPEVAL  ( NN  , X  ,  Y  , FDP  , XX  , F  ) 


AK2 ( I  )  =F 
280  CONTINUE 


CALL  STRESS  (ND  ,N  ,CHVD ,DT ,ESTI , FR , S I GMAO , S IGMAV , S IGZER , AKO , RA ( N58 ) 


DO  290  1  =  1  , ND 

GI (  I  )  =  1000 . *AK2 ( I ) *DSQRT ( SIGMAO ( I ) ) 
290  CONTINUE 

T =DT*KKK 
KKK=KKK+1 

IF(T.LT.TMAX )  GO  TO  180 
STOP 


300  FORMAT  ( / /5X , ' TOTAL  STORAGE  USED  FOR  THIS  PROBLEM 
310  FORMAT  ( 1H1 , ///10X ,' STORAGE  EXCEEDED',//) 

320  FORMAT  (8F10.0) 

330  FORMAT  (8F10.4) 

340  FORMAT  (10F5.0) 

350  FORMAT  ( 3X , F7 . 4 , F7 . 3 ,8F7 . 3  ) 

360  FORMAT  (10F5.0) 


FORMAT  (10F5.1) 

FORMAT  (3F10.0) 

FORMAT  (3F10.4) 

FORMAT  (FIO.O) 

FORMAT  (F10.4) 

FORMAT  (4F10.0) 

FORMAT  (4F10.4) 

FORMAT  (FIO.O, 2F10.0) 

FORMAT  ( FIO . 4 , 2F 10 . 4 ) 

FORMAT  (4F10.0) 

FORMAT  (4F10.4) 

FORMAT  (2F10.0) 

FORMAT  (2F10.4) 

FORMAT  (4F10.0) 

FORMAT  (4F10.4) 

FORMAT  (110) 

FORMAT  (110) 

FORMAT  (1H1.///15X, 'SOIL  PROPERTIES  ARE',///) 

FORMAT  (5X, 'TOTAL  DEPTH  OF  SAND  STRATA  » '  .F15.4/5X,  'ACCELE 

1RATI0N  DUE  TO  GRAVITY  =' ,F15 .4/5X  ,' DAMPING  COEFFICIENT 

2  =' .F15.4/5X, 'UNIT  WEIGHT  OF  SATURATED  SAND  =',F15. 

34 /5X, 'UNIT  WEIGHT  OF  WATER  = '  ,  FI  5 . 4/5X  ,  ' SPEC  I F I C  GRA 

4VITY  OF  SAND  =' ,F15 .4/5X  ,' ANGLE  OF  FRICTION 

5  =' .F15.4/5X, 'COEFFICIENT  OF  EARTH  PRESSURE  = '  ,F15 .4/5X  ,  ' 

6V0ID  RATIO  ='  ,F15.4/5X, 'UNIT  WEIGHT  OF  DRY 

7  SAND  =' .F15.4/5X, 'HEIGHT  OF  EACH  LAYER 

8'  .F15.4/5X,  'RELATIVE  DENSITY  OF  SAND  = '  ,F15 .4/// ) 

FORMAT  ( 1H1 ) 

FORMAT  ( //  ,  1 5X  ,  ' AT  TIME  T  =  0.0  SEC ' //3X  ,' TOTAL  STRESS  '  ,3X  ,' EFFECTIV 
1ESTRESS' ,3X, 'PORE  PRESSURE',//) 

FORMAT  ( 3 ( 2X  ,  F 1 5 . 4  )  ) 

FORMAT  (1H1.21X, 'CYCLE  NUMBER  IS*  ,15, /24X  ,  ' . ’,/) 

FORMAT  (10X, 'DAMPING  RATIO  =  '  ,F10 . 5/ ,  10X , ' LAYER  FIRST  TO  LIQUEFY 
1=  ',15,/) 

FORMAT  (2X,  'TIME' , 8X , ' SHE AR ( 2 ) ' , 2X , ' SHE AR ( 3 ) ' , 2X , ’ SHEAR ( 4 ) ' ,2X, ' SH 
1EAR ( 5 ) ' , 2 X , 'SHEAR (6) ' ,2X, ' SHEAR ( 7 ) ' ,2X, ' SHEAR (8) ’ , 2X , ’ SHE AR ( 9 ) ' ,1X 
2,  'SHEAR(IO) ' ,1X, ' SUR .  ACCE’,/) 


o  o  o  o 
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SUBROUTINE  SPLINE  ( NN , X , Y , FDP , A , B , C  , R  ) 

IMPLICIT  REAL*8 ( A-H  ,0-Z  ) 

C 

COMMON  RA { 18000  ) 

DIMENSION  X ( NN  )  ,  Y ( NN  )  ,  A ( NN  )  ,  B  (  NN  )  ,  C(NN),  R ( NN  )  ,  FDP(NN) 

THIS  IS  A  PART  OF  SPLINE  INTERPOLATION  USED  TO  CALCULATE 
THE  VALUE  OF  K2 

THE  SPLINE  HAS  PARABOLIC  RUNOUTS 

ALAMDA  =  1 .0 
NM2=NN-2 
NMl =NN - 1 
C 

C(U=X(2)-X(1) 

C 

DO  10  1=2, NMl 
C(I)=X(I+1)-X(I) 

A ( I ) =C (  I - 1 ) 

B  ( I  ) =2 . * ( A ( I )+C ( I ) ) 

R( I  )  =6 . * ( ( Y(I  +  1  ) - Y ( I ) )/C( I ) - ( Y ( I )-Y( 1-1 ) )/C( 1-1 )) 

10  CONTINUE 

C 

B(2)=B(2)+ALAMDA*C(1) 

B(NM1)=B(NM1)+ALAMDA*C(NM1  ) 

C 

DO  20  1=3, NMl 
T  =  A( I  )/B(  I  —  1 ) 

B ( 1 )  =  B ( I ) -T*C ( I - 1 ) 

R( I )=R( I  )-T*R(  1-1 ) 

20  CONTINUE 
C 

FDP (NMl )  =  R(NM1 )/B(NMl ) 

C 

DO  30  1=2, NM2 
MMI-NN-I 

FDP (NMl  )  =  ( R { NMl ) -C ( NMl )*FOP (NMI  +  1 ) ) /B { NMl ) 

30  CONTINUE 
C 

FDP ( 1  )  =AL AMDA*FDP ( 2 ) 

FDP (NN  )=ALAMDA*FDP (NMl ) 

C 

RETURN 

END 


SUBROUTINE  TLOAD  ( NTP ,FF , NT ,ND , SM, I  COUNT , ACC  ) 
IMPLICIT  REAL*8 ( A-H  ,0-Z ) 

COMMON  RA ( 18000  ) 

COMMON  /THREE/  OELT  ,TMAX ,AA  .ANOMEG 
DIMENSION  FF ( ND  ,NTP  ) ,  SM(ND.ND),  ACC(NT) 

SUBROUTINE  CALCULATES  LOAD  HISTORY 

I F ( I  COUNT  .  NE  .  1 )  GO  TO  40 

NT1=NT+1 

PI=4.*DATAN(1.D0) 

OMEGA=ANOMEG*P I 

CALCULATE  LOAD  AT  EACH  STEP  OF  DELT  =  0.01  SEC 


J  =  1 

DO  10  1=1, ND 
FF ( I , J  )  =0 . 0 

DO  30  1=1, ND 
DO  20  J  =  1  ,NT 
AJ  =  J 

T  =  AJ*DELT 

FF(I,J  +  1)=-SM(I,I )*AA*OSIN( OMEGA*T ) 

CONTINUE 

CONTINUE 

GO  TO  80 
CONTINUE 

READ  (5,90)  ( ACC ( I ) , I  =  1 , NT  ) 

CALCULATE  LOAD  AT  EACH  STEP  OF  DELT  =  0.01  SEC 


J  =  1 

DO  50  1  =  1, ND 
FF ( I , J  )=0.0 

DO  70  1  =  1, ND 
DO  60  J=1,NT 

FF(I,J+1)=-SM(I,I)*ACC(J) 

CONTINUE 

CONTINUE 

CONTINUE 


RETURN 


90  FORMAT  (8F10.0) 


SUBROUTINE  STIFF  (EKI,ND,SK) 


IMPLICIT  REAL*8(A-H,0-Z) 
DIMENSION  EKI(ND),  SK  (ND  ,ND  ) 

DO  10  1=1, ND 
DO  10  J  =  1  ,ND 
SK ( I , J ) =0 . 0 

SK ( 1 ,1 )=EKI (1  ) 

DO  20  1=2, ND 

SK(I,I)=EKI(I-1)+EKI(I) 

CONTINUE 

DO  30  J  =  2  ,  ND 
I  =J  - 1 

SK ( I ,J  )=-EKI  (  I  ) 

CONTINUE 

NDM1=ND-1 
DO  40  J=1 ,NDM1 
I  =J  +  1 

SK ( I ,J )=  —  EKI  (J  ) 

CONTINUE 


RETURN 


o  o  o 


c 


r.'* 

V 


SUBROUTINE  DAMP  ( WN ,ZETA , ND  ,  SK  ,  SC  ) 

IMPLICIT  REAL*8 ( A-H  ,0-Z  ) 

DIMENSION  SK(ND.ND),  SC(ND,ND) 

CALCULATE  CONSTANT  A1  IN  RAYLEIGH  DAMPING 


AI =2 . *ZETA/WN 
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SUBROUTINE  D YNA ( B IG , SM , $K , SC , FF  ,G I  ,H, 

UDI ,UVI ,UAI ,UD,UV,UA,SHEAR,F,SST,A,B,X) 


IMPLICIT  RE AL*8 ( A-H  , O-Z  ) 

COMMON  RA ( 18000  ) 

COMMON  /ONE/  KKK 

COMMON  /TWO/  I L I  0 

COMMON  /AAA/  ND ,LAYDRY ,NN  ,N  ,NM 

COMMON  /BBB /  NT , ND I V  ,NTP  ,  ND  I  VP 

COMMON  /FOUR/  TSTEP , DELTA , ALPHA ,TF I N 

DIMENSION  SM(ND ,ND ) ,  SK(ND.ND),  SC(ND,ND),  A(ND  ,ND ) ,  B ( ND  ) ,  X(ND), 
1  UOI(ND),  UVI(ND),  UAI(ND),  UD(ND),  UV(ND),  UA(ND),  SST(ND),  GI(ND 
2),  BIG (ND ) 

DIMENSION  H ( ND ) ,  SHEAR ( ND ) ,  F(ND.NDIVP),  FF ( ND  ,NTP  ) 

THIS  SUBROUTINE  SOLVES  THE  DYNAMIC  EQUATION  OF  MOTION  TO  GET 
THE  SHEAR  STRAINS 

NT1=NDI V+l 

CORELATE  THE  LOADS  WITH  TOTAL  LOAD  GIVEN  BY  FF 

DO  10  1  =  1, ND 
DO  10  J  =  1  ,NT1 
F ( I , J  )  =0 . 0 
DO  30  1  =  1, ND 
IF(I.LE.ILIQ)  GO  TO  30 
DO  20  J  =  1  ,NT1 
L=J+50*( KKK-1  ) 

F(I,J)=FF(I,L) 

CONTINUE 

CONTINUE 

ONLY  FOR  THE  FIRST  CYCLE  GO  THROUGH  THE  PROCEDURE 
OF  EVALUATING  THE  INITIAL  ACCELERATION  AND  VARIOUS 
CONSTANTS  IN  NEWMARK'S  SCHEME  FOR  STEP-BY-STEP 
INTEGRATION. 

IF (KKK . NE . 1  )  GO  TO  80 

DO  40  1  =  1, ND 
DO  40  J  =  1  ,ND 
A ( I , J  )  =SM ( I ,  J  ) 


DO  50  1=1, ND 
B { I ) =F ( I , 1 ) 


115 


DO  60  1  =  1, NO 
DO  60  J  =  1 ,  ND 

60  B ( I ) =B ( I ) -SC ( I ,J ) *  U  V  I (J  ) -SK ( I  ,J  )  *UD  I  (J) 

C 

CALL  TRI  (A,ND  ) 

CALL  REDUCE  (A,X,B,ND) 

C 

C 

DT  =  TSTEP 

A0=1./(ALPHA*DT*DT) 

A1=DELTA/(ALPHA*DT) 

A2  =  l . /( ALPHA*DT  ) 

A3=1./(2.*ALPHA)-1. 

A4=DELTA/ ALPHA-1. 

A5  =  DT*0. 5* (DELTA/ALPHA-2.  ) 

A6=DT*(1. -DELTA) 

A7=DELTA*DT 

C 

DO  70  1  =  1, ND 

70  UAI ( I  )  =  X  ( I  ) 

C 

80  CONTINUE 
C 

DO  150  L  =  1  ,NDIV 
DO  90  1  =  1, ND 
DO  90  J  =  1  ,ND 

90  A(I,J)=SK(I,J)+A0*SM(I,J)+A1*SC(I,J) 

DO  110  1  =  1, ND 
B  ( I )  =F ( I , L  +  l ) 

DO  100  J=1,ND 

100  B ( I ) =B ( I )+SM ( I , J )*(AO*UDI (J )+A2*UVI (J )+A3*UAI (J ) )+SC( I ,J )*(A1*UDI ( 
1J )+A4*UVI (J )+A5*UAI (J ) ) 

110  CONTINUE 
C 

CALL  TRI  (A.ND) 

CALL  REDUCE  (A,X,B,ND) 

C 

AL  =  L 
T=DT*AL 
DO  120  1=1, ND 
UD ( I ) =X  ( I  ) 

C 

U A ( I )=A0*( UD ( I  ) - UD I ( I ) )-A2*UV I ( I )-A3*UAI  ( I ) 

C 

UV ( I  )  =  U V  I ( I  )  +A6*UAI ( I  )  +  A7*UA( I ) 

C 

120  CONTINUE 

C  CALCULATE  SHEAR  STRAINS  FOR  EACH  LAYER 
C 

ND 1  =  ND- 1 
DO  130  1  =  1, ND1 


$ST(I)=(UD(I)-UD(I+1))/H(I) 

SHEAR ( I ) =GI ( I )*SST { I ) 

SST( I  )=DABS(SST( I )) 

130  CONTINUE 

SST(ND)=UD(ND)/H(ND) 

SHEARING  )=GI (ND)*SST (ND) 

SST(ND)=DABS(SST(ND)) 

CALCULATE  THE  MAXIMUM  SHEAR  STRAIN  FOR  A  PAERIOD 
OF  DT  =  0.5  SEC 

CALL  AMAX  ( ND , L , SST  , B  IG  ) 


DO  140  1=1, ND 
UDKI)-UD(I) 

UAI ( I ) =UA (  I  ) 

UVKI)-UV(I) 

140  CONTINUE 

I C  =  L  +  50* ( KKK - 1  ) 

A I C  =  I  C 
TE=DT*AIC 

WRITE  (6,160)  TE,(SHEAR(IJ),IJ=2,10),UAI(1) 
150  CONTINUE 
RETURN 


160  FORMAT  ( F10 . 5 , 9F10 . 3 , F10 . 5  ) 


SUBROUTINE  AMAX  (ND ,L ,SST  ,BIG) 


IMPLICIT  REAL*8 { A-H  ,0-Z  ) 

DIMENSION  SST(ND),  BIG(NO) 
IF(L.NE.l)  GO  TO  20 

DO  10  1  =  1, ND 
BIG ( I )  =0 .0 
CONTINUE 

CONTINUE 

DO  30  I  =1 ,  ND 

IF ( SST ( I  ) . GT  . B  I G  ( I ) )  BIG ( I  )  =  SST ( I  ) 
CONTINUE 


RETURN 


OOO  OOO  OOOOO 
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SUBROUTINE  SPEVAL  ( NN , X , Y , FDP , XX  , F  ) 


IMPLICIT  RE AL *8 ( A-H , 0-Z  ) 

THIS  SUBROUTINE  EVALUATES  THE  CUBIC  SPLINE  INTERPOLATED 
VALUE  OF  K2 

F=INTERPOLATED  VALUE 
DIMENSION  X ( NN  )  ,  Y  ( NN  )  ,  FDP(NN) 

FIND  THE  PROPER  INTERVAL 

NM1=NN-1 
DO  10  1=1, NM1 
I F ( XX , LE . X ( I  +  1 ) )  GO  TO  20 
10  CONTINUE 
RETURN 

EVALUATE  THE  CUBIC 

20  DXM=XX-X {  I  ) 

DXP  =  X ( 1  +  1  )  -XX 
DEL  =  X ( 1  +  1  )-X  (  I  ) 

C 

F  =  FDP  (  I  )*DXP*(DXP*DXP/DEL-DEL)/6.+FDP(I+l)*DXM*(DXM*DXM/DEL-DEL)/6 
1 . +Y ( I )*DXP/DEL+Y( 1+1 )*DXM/DEL 
RETURN 
END 


IMPLICIT  RE AL *8 ( A-H  , 0-Z  ) 

COMMON  RA ( 18000  ) 

COMMON  /FIVE/  C1,C2,C3,C4 

DIMENSION  CHVD(ND),  EPSVD(ND),  DELVD(ND),  BIG (NO ) 

THIS  SUBROUTINE  CALCULATES  VALUE  OF  THE  SOURCE  TERM  IN  THE 
DIFFUSION  EQUATION 

DO  10  1=1, ND 

DELVDU  )=C1*(BIG(  I  )-C2*EPSVD(I )  )+C3*{  EPSVD  ( I  )**2  )/(  BIG  ( I  )+C4*EPSVD 
(I)) 

CHVD ( I ) =DELVD (  I  ) 

CONTINUE 

DO  20  1=1, ND 

EPSVD ( I )=DELVD(  I  )+EPSVD(  I  ) 


RETURN 


OOOOC1  o  <->  ooooo 
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SUBROUTINE  DIFF  ( I FL AG , N ,FR ,Q , BETA  , NM , A  , B  ,C  , R  ,T  , BB  ) 

C 

IMPLICIT  REAL*8 ( A-H  ,0-Z  ) 

C 

DIMENSION  FR { N  ) ,  T  {  N  ) ,  Q  (  N  ) ,  BETA ( N  ) ,  A ( N  ) ,  B ( N  )  ,  C(N),  R  (  N  ) ,  BB ( N 
1) 

C 

NM1=N-1 

IF ( IFLAG.NE . 1 )  GO  TO  30 

TEMPORARY  ARRAY  IS  DENOTED  BY  T(l) 

THIS  IS  EXPLICIT  METHOD  SCHEME 


DO  10  1=3, NM1 

T {  I  )  =FR (  I  )  +  BETA( I )*(FR( 1-1 )-2.*FR( I )  +  FR ( 1  +  1 ) )+Q ( I  ) 
10  CONTINUE 

T(2)=FR(2)+BETA(2)*(-2.*FR(2)+2.*FR(3))+Q(2) 

T(N)=FR(N)+BETA{N)*(-2.*FR(N)+FR(NM1))+Q(N) 

DO  20  1=2, N 
FR ( I ) =T ( I  ) 

20  CONTINUE 
RETURN 

30  CONTINUE 

THIS  IS  CRANK-NICOLSON  IMPLICIT  SCHEME 


DO  40  1=2, N 
A ( I ) =- BETA (  I  ) 

B(I)=2.*(1.+BETA(I)) 

C ( I ) =-BETA ( I  ) 

40  CONTINUE 
C  START  THE  MAIN  LOOP 
C 

DO  50  1=3, NM1 

R(  I )  =  B E T A (  I )*{FR{ I-1)  +  FR{ 1  +  1) ) +  2 . * ( 1 . -BETA ( I ))*FR(I  )+Q(I  ) 
50  CONTINUE 
C 

R(2)=BETA(2)*FR(3)*2.+2.*(1.-BETA(2))*FR(2)+Q(2) 
R(N)=BETA(N)*FR(N-1)+2.*(1.-BETA(N) )*FR(N)+Q(N) 

C 

CALL  SOLVE  ( N  ,  A  ,B  ,  C  , FR  ,  R  ,  BB  ) 

RETURN 

END 


SUBROUTINE  STRESS  ( ND ,N ,CHVD ,DT , EST I , FR  ,  S IGMAO , S IGMAV , SIGZER ,AKO ,B 
IB) 

C 

IMPLICIT  RE  AL*8 ( A-H  ,  O-Z  ) 

C 

COMMON  RA ( 18000  ) 

COMMON  /ONE/  KKK 
COMMON  /TWO/  I L  IQ 
C 

DIMENSION  BB(ND),  S IGMAV ( ND  )  ,  S IGMAO ( ND  ) ,  ESTI(ND),  CHVD(ND),  SIGZ 
1ER(ND),  FR ( N  ) 

C 

C 

C  DEFINE  INCREASE  IN  THE  PORE  PRESSURE  AT  THE  CENTER  OF 
C  EACH  LAYER 

C 

TT  =KKK*DT 
C 

C  BB=INCREASE  IN  THE  PORE  PRESSURE 

C 

BB ( 1 ) =0 . 0 
J  =  N 

DO  10  I =2  ,ND 
BB ( I ) =FR ( J  ) 

J  =  J-2 

10  CONTINUE 
C 
C 

DO  20  I  =  1  ,  ND 

S IGMAV ( I )=SIGZER(  I  )  — BB(  I  ) 

IF(SIGMAV(I).LE.O.O)  I L I Q= I 
20  CONTINUE 
C 
C 

DO  50  1  =  1, ND 
CHECK=1.E-5*SIGZER(  I  ) 

IF ( SIGMAV { I  )  -  CHECK  )  30,40,40 
30  SIGMAV ( I  )=CHECK 
40  CONTINUE 
50  CONTINUE 
C 

WRITE  (6,80) 

C 

WRITE  (6,90)  TT,(SIGMAV(I),BB(I),I=1,ND) 

C 

DO  60  1  =  1, ND 

SIGMAO ( I ) = { SIGMAV ( I ) +2 . *AKO*S IGMAV ( I ) ) / 3 . 

60  CONTINUE 
C 

C  CONVERT  SIGMAV  INTO  ESTI  FOR  THE  NEXT  CYCLE 
C 


DO  70  1  =  1,  ND 
ESTI ( I ) =SIGMAV ( I  ) 
CONTINUE 


RETURN 

FORMAT  ( ///25X , 1  NEW  EFFECTIVE  STRESSES  AND  PORE  PRESSURES  ARE',/) 
FORMAT  ( / / 3 0 X , ' TIME= 1 , F10 . 4//1 { 25X , F10 . 4 , 5X  , F 10 . 4  ) ) 


SUBROUTINE  TRI  ( A , ND  ) 

IMPLICIT  REAL*8{ A-H  ,0-Z  ) 

DIMENSION  A(ND  ,ND  ) 

ND 1 =ND - 1 
TINY=1.E-15 

DO  30  1=1, ND1 

IF (DABS(A( I ,1  ) ) .LT.TINY  )  GO  TO  40 
J1=I+1 


DO  20  J  =  J1  ,ND 
I F ( A ( J , I  ) . EQ  .  0  .  )  GO  TO  20 
FAC=A ( J , I ) /A (  I  ,  I  ) 

DO  10  K  =  J1  ,ND 

A(J,K)=A(J,K)-A(I,K)*FAC 

CONTINUE 

CONTINUE 

RETURN 

WRITE  (6,50)  I.Afl.I) 

STOP 


FORMAT  (//5X, 'SMALL  PIVOT ' , 15  ,F10 .4  , /) 
END 
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SUBROUTINE  REDUCE  (A,X,B,ND) 
c  ================================ 

IMPLICIT  REAL*8 ( A-H  ,0-Z  ) 

C 

DIMENSION  A(ND'ND),  B ( NO  ) ,  X(ND) 
C 

ND1=ND-1 

C 

DO  20  1  =  1, ND1 
J 1 = I +1 
C 

no  10  J  =  J1  ,ND 

10  B(J)=B(J)-B(I)*A(J,I)/A(I,I) 

20  CONTINUE 
C 

X(ND)=B(ND)/A(ND,ND) 

C 

DO  40  1  =  1, ND1 
I  B  =  ND - 1 
J 1 = I B+l 
DO  30  J= J 1 , ND 

30  B(IB)=B(IB)-A(IB,J)*X(J) 

X( IB)  =  B( IB)/A( IB,  IB) 

40  CONTINUE 
RETURN 
END 


ooo  o  ooooooooo 
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SUBROUTINE  SOLVE  ( N , A  , B , C  , X  , Q  , BB  ) 

IMPLICIT  REAL*8 { A-H  ,0-Z  ) 

C 

DIMENSION  A ( M ) ,  B(N),  C(N),  X(N),  Q(N),  BB  (  N  ) 

THIS  SUBROUTINE  SOLVES  TRIDIAGONAL  SYSTEM  UDSING 
GAUSS  ELIMINATION. 

( N  —  1 )  IS  THE  TOTAL  NUMBER  OF  EQUATIONS  TO  BE  SOLVED. 
BB  IS  DUMMY  ARRAY. 

FORWARD  ELIMINATION. 

DO  10  1=2, N 
BB ( I  )  =  B ( I ) 

10  CONTINUE 

DO  20  1=3, N 
T=A(I)/BB(I-1) 

BB(I)=BB(I)-C(I-1)*T 
Q  ( I ) =Q ( I ) -Q ( I  - 1 ) *T 
20  CONTINUE 

BACK  SUSTITUTE 

X(N)=Q(N)/BB(N) 

C 

NM2=N-2 
DO  30  1=1, NM2 
J  =  N- 1 

X(J )  =  (Q(J )-C(J  ) *X ( J  +  1 )  )/BB  (J  ) 

30  CONTINUE 
C 

RETURN 

END 


Appendix  E 
SAMPLE  OUTPUT 


The  sample  output  for  EXAMPLE-1  using  explicit  method 
is  given  below. 

THIS  IS  THE  DYNAMIC  RESPONSE  ANALYSIS  OF  LAYERED  SAND 
STRATUM. 


TOTAL  NUMBER  OF  LAYERS  =  10 

THIS  IS  THE  DYNAMIC  RESPONSE  ANALYSIS  OF  LAYERED  SAND 
STRATUM. 

TOTAL  NUMBER  OF  LAYERS  =  10 

NUMBER  OF  DRY  LAYERS  =  1 

NUMBER  OF  DATA  POINTS  FOR  INTERPOLATION  =  10 

NUMBER  OF  DIVISIONS  OF  EACH  LAYER  =  2 

TOTAL  NUMBER  OF  SPATIAL  DIVISIONS  =  18 

NUMBER  OF  DIVISIONS  OF  EXCITATION  PERIOD  =  1000 

NUMBER  OF  DIVISIONS  OF  EACH  CYCLE  PERIOD  =  50 


TOTAL  STORAGE  USED  FOR  THIS  PROBLEM  =  12511 

THIS  IS  INPUT  DATA  FOR  THIS  PROBLEM. 


50.0000  32.2000  0.0200  122.0000  62.4000  2.6500 

40.000  0.4500 

0.0001  0.001  0.010  0.100  0.200  0.400  0.500  0.700 

0.800  1.00 


0.4300 

0.6200 

0.0025 

0.0000 

0.0003 

45.0000 

0.5000 

10.0000 

4.0000 

0.0100 

0.8050 

0.0100 

0.5000 

0.2500 

0.5000 

0.0000 

0  .0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.8000 

0.7900 

0.4500 

0.6300 

1 


SOIL  PROPERTIES  ARE 


TOTAL  DEPTH  OF  SAND  STRATA 
ACCELERATION  DUE  TO  GRAVITY 
DAMPING  COEFFICIENT 
UNIT  WEIGHT  OF  SATURATED  SAND 
UNIT  WEIGHT  OF  WATER 
SPECIFIC  GRAVITY  OF  SAND 
ANGLE  OF  FRICTION 
COEFFICIENT  OF  EARTH  PRESSURE 
VOID  RATIO 

UNIT  WEIGHT  OF  DRY  SAND 
HEIGHT  OF  EACH  LAYER 
RELATIVE  DENSITY  OF  SAND 


50.0000  ft. 
32.2000  ft/sec 
0.0200 

122.0000  lb/ft 
62.4000  lb/ft 
2.6500 

40.0000  deg. 
0.3572 
0.7275 

95.7212  lb/ft 
5.0000  ft. 
0.4500 


AT  TIME  T=0 . 0  SEC 


TOTAL  STRESS 


EFFECTIVE  STRESS  PORE  PRESSURE 


239.3030  239 

783.6061  627 

1393.6061  925 

2003.6061  1223 

2613.6061  1521 

3223.6061  1819 

3833.6061  2117 

4443.6061  2415 

5053.6061  2713 

5663.6061  3011 


3030 

0.0000 

6061 

156.0000 

6061 

468.0000 

6061 

780.0000 

6061 

1092.0000 

6061 

1404.0000 

6061 

1716.0000 

6061 

2028.0000 

6061 

2340.0000 

6061 

2652.0000 

CYCLE  NUMBER  IS 


1 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  =  0 

TIME  SHEAR ( 2 )  SHEAR ( 3 )  SHEAR { 4 )  SHEAR  ( 5 


SHEAR ( 6 ) 

0.01000  -0 

0.02000  -0 

0.03000  -0 

0.04000  -0 

0.05000  -0 

0.06000  -0 

0.07000  -0 

0.08000  -1 

0.09000  -4 

0.10000  -7 

0.11000  -12 

0.12000  -17 

0.13000  -22 

0.14000  -27 

0.15000  -31 

0.16000  -34 

0.17000  -36 

0.18000  -38 

0.19000  -38 

0.20000  -38 

0.21000  -38 

0.22000  -36 

0.23000  -33 

0.24000  -30 

0.25000  -25 

0.26000  -19 

0.27000  -12 

0.28000  -4 

0.29000  3 

0.30000  12 

0.31000  21 

0.32000  29 

0.33000  37 

0.34000  43 

0.35000  48 

0.36000  52 

0.37000  55 

0.38000  57 

0.39000  57 

0.40000  56 

0.41000  54 

0.42000  50 

0.43000  45 


.000  -0.000 

.000  -0.001 

.002  -0.009 

.014  -0.055 

.069  -0.237 

.260  -0.793 

.787  -2.137 

.969  -4.791 

.172  -9.175 

.625  -15.355 

.229  -22.962 

.511  -31.335 

.802  -39.770 

.527  -47.669 

.403  -54.570 

.407  -60.151 

.631  -64.269 

.144  -66.944 

.947  -68.253 

.989  -68.209 

.196  -66.725 

.499  -63.664 

.849  -58.894 

.199  -52.315 

.498  -43.861 

.698  -33.524 

.811  -21.398 

.956  -7.735 

.613  7.031 

.531  22.302 

.370  37.401 

.698  51.674 

.156  64.580 

.510  75.736 

.664  84.917 

.618  92.021 

.403  97.006 

.024  99.840 

.439  100.468 

.577  98.821 

.372  94.839 

.795  88.501 

.874  79.863 


-0.000  -0 
-0.004  -0 

-0.033  -0 

-0.175  -0 

-0.665  -1 

-1.955  -4 

-4.653  -8 

-9.287  -15 

-16.030  -24 

-24.627  -35 

-34.569  -46 

-45.306  -59 

-56.288  -72 

-66.901  -85 

-76.481  -97 

-84.461  -107 

-90.492  -114 

-94.440  -120 

-96.300  -122 

-96.089  -122 

-93.769  -118 

-89.207  -112 

-82.218  -103 

-72.637  -90 

-60.400  -74 

-45.581  -55 

-28.416  -33 

-9.329  -9 

11.070  15 

32.027  41 

52.727  67 

72.390  91 

90.347  114 

106.072  134 

119.179  151 

129.401  164 

136.565  173 

140.565  178 

141.335  179 

138.833  176 

133.037  168 

123.970  156 

111.741  141 


.001  -0.004 
.015  -0.051 
.107  -0.315 
.492  -1.241 
.630  -3.547 
.189  -7.934 
.781  -14.694 
.639  -23.626 
.554  -34.345 
.106  -46.546 
.907  -59.972 
.580  -74.245 
.620  -88.840 
.362  -103.153 
.090  -116.514 
.139  -128.196 
.964  -137.499 
.189  -143.849 
.634  -146.855 
.248  -146.288 
.991  -142.030 
.761  -133.999 
.389  -122.106 
.727  -106.298 
.741  -86.654 
.605  -63.474 
.721  -37.310 
.697  -8.908 
.709  20.876 

.658  51.160 
.282  81.073 
.750  109.762 
.299  136.387 
.257  160.139 
.056  180.274 
.245  196.168 
.496  207.342 
.609  213.488 
.484  214.449 
.101  210.202 
.494  200.842 
.763  186.575 
.113  167.737 


0.44000 

39.704 

69.082 

96.589 

121  .893 

144.809 

0.45000 

32.438 

56.431 

78.908 

99.617 

118.426 

0.46000 

24.292 

42.297 

59.239 

74.957 

89.362 

0.47000 

15.540 

27.162 

38.242 

48.704 

58.495 

0.48000 

6.506 

11.569 

16.639 

21.709 

26.753 

0.49000 

-2.467 

-3.921 

-4.843 

-5.182 

-4.939 

0.50000 

-11.048 

-18.780 

-25.526 

-31.176 

-35.706 

SHEAR ( 7  ) 

SHEAR { 8 ) 

SHEAR ( 9 ) 

SHEAR ( 10 ) 

SUR.  ACCE 

0.015 

-0.054 

-0.186 

-0.621 

-0.10089 

0.162 

-0.478 

-1.301 

-3.202 

-0.20019 

0.842 

-2.043 

-4.439 

-8.424 

-0.29629 

2.810 

-5.683 

-10.170 

-16.072 

-0.38742 

6.872 

-11.855 

-18.309 

-25.889 

-0.47111 

3.393 

-20.372 

-28.509 

-37.712 

-0.54276 

2.179 

-30.864 

-40.589 

-51.305 

-0.59340 

2.852 

-43.110 

-54.344 

-66.373 

-0.60794 

5.155 

-56.915 

-69.463 

-82.616 

-0.56708 

8.879 

-71.970 

-85.619 

-99.709 

-0.45545 

3.702 

-87.915 

-102.487 

-117.295 

-0.27357 

9.212 

-104.387 

-119.691 

-134.992 

-0.04477 

4.953 

-120.965 

-136.804 

-152.389 

0.19206 

0.395 

-137.133 

-153.351 

-169.017 

0.40246 

4.887 

-152.287 

-168.769 

-184.331 

0.57413 

7.694 

-165.729 

-182.380 

-197.678 

0.71713 

8.053 

-176.667 

-193.394 

-208.279 

0.84705 

5.239 

-184.272 

-200.935 

-215.241 

0.96818 

8.639 

-187.766 

-204.113 

-217.628 

1.07188 

7.815 

-186.509 

-202.139 

-214.584 

1.14754 

2.510 

-180.079 

-194.460 

-205.492 

1.19142 

2.616 

-168.316 

-180.862 

-190.113 

1.20482 

8.139 

-151.312 

-161.501 

-168.648 

1.18694 

9.198 

-129.380 

-136.857 

-141 .684 

1.13221 

6.063 

-103.016 

-107.633 

-110.055 

1.03369 

9.192 

-72.858 

-74.629 

-74.666 

0.88732 

9.253 

-39.662 

-38.664 

-36.387 

0.69320 

7.063 

-4.254 

-0.554 

3.972 

0.45546 

6.489 

32.497 

38.875 

45.602 

0.18298 

0.513 

69.710 

78.752 

87.627 

-0.10974 

4.133 

106.487 

118.140 

129.074 

-0.40396 

6.467 

141.894 

156.044 

168.896 

-0.68070 

6.605 

174.954 

191.428 

206.014 

-0.92528 

3.625 

204.667 

223.244 

239.352 

-1.13000 

6.646 

230.061 

250.470 

267.862 

-1.29394 

4.896 

250.257 

272.154 

290.552 

-1.42030 

7 .764 

264.525 

287.477 

306.550 

-1.51261 

4.822 

272.328 

295.814 

315.180 

-1.57172 

5.841 

273.359 

296.795 

316.039 

-1.59488 

0.782 

267.562 

290.352 

309.054 

-1.57727 

9.793 

255.133 

276.719 

294.481 

-1.51462 

3.199 

236.500 

256.402 

272.874 

-1.40520 

134 


CYCLE  NUMBER  IS  2 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  =  0 

TIME  SHE AR ( 2  )  SHEAR ( 3  )  SHEAR ( 4  )  SHEAR ( 5  ) 

SHEAR ( 6 ) 


0.51000 

-12.885 

-21.457 

0.52000 

-17.699 

-29.855 

0.53000 

-21.782 

-37.396 

0.54000 

-25.141 

-43.918 

0.55000 

-28.014 

-49.388 

0.56000 

-30.626 

-53.958 

0.57000 

-32.992 

-57.861 

0.58000 

-34.973 

-61.248 

0.59000 

-36.489 

-64.133 

0.60000 

-37.625 

-66.470 

0.61000 

-38.534 

-68.260 

0.62000 

-39.293 

-69.546 

0.63000 

-39.849 

-70.317 

0.64000 

-40.075 

-70.458 

0.65000 

-39.837 

-69.789 

0.66000 

-39.030 

-68.155 

0.67000 

-37.580 

-65.474 

0.68000 

-35.454 

-61.732 

0.69000 

-32.678 

-56.960 

0.70000 

-29.333 

-51.244 

0.71000 

-25.535 

-44.729 

0.72000 

-21.408 

-37.610 

0.73000 

-17.066 

-30.090 

0.74000 

-12.607 

-22.330 

0.75000 

-8.106 

-14.433 

0.76000 

-3.605 

-6.444 

0.77000 

0.897 

1.628 

0.78000 

5.432 

9.791 

0.79000 

10.038 

18.043 

0.80000 

14.729 

26.365 

0.81000 

19.484 

34.710 

0.82000 

24.241 

42.995 

0.83000 

28.911 

51.100 

0.84000 

33.395 

58.871 

0.85000 

37.586 

66.120 

0.86000 

41.375 

72.642 

0.87000 

44.647 

78.231 

0.88000 

47.287 

82.704 

0.89000 

49.188 

85.908 

0.90000 

50.274 

87.737 

0.91000 

50.503 

88.131 

0.92000 

49.876 

87.080 

-29.131 

-35.794 

-41.451 

-40.833 

-50.575 

-59.057 

-51.537 

-64.234 

-75.456 

-61.080 

-76.593 

-90.456 

-69.330 

-87.516 

-103.937 

-76.260 

-96.934 

-115.848 

-82.010 

-104.889 

-126.177 

-86.839 

-111.538 

-134.912 

-90.976 

-117.076 

-142.049 

-94.480 

-121.615 

-147.605 

-97.246 

-125.103 

-151.598 

-99.108 

-127.364 

-154.002 

-99.953 

-128.217 

-154.717 

-99.721 

-127.561 

-153.612 

-98.352 

-125.371 

-150.610 

-95.748 

-121.639 

-145.721 

-91.807 

-116.339 

-139.020 

-86.481 

-109.456 

-130.596 

-79.811 

-101.029 

-120.530 

-71.915 

-91.171 

-108.910 

-62.968 

-80.055 

-95.848 

-53.174 

-67.863 

-81.466 

-42.741 

-54.771 

-65.899 

-31.850 

-40.935 

-49.294 

-20.632 

-26.507 

-31.827 

-9.168 

-11.632 

-13.706 

2.488 

3.550 

4.850 

14.287 

18.912 

23.638 

26.169 

34.336 

42.475 

38.065 

49.708 

61.182 

49.892 

64.896 

79.559 

61.544 

79.741 

97.371 

72.869 

94.036 

114.352 

83.661 

107.528 

130.215 

93.666 

119.933 

144.666 

102.613 

130.951 

157.410 

110.236 

140.295 

168.165 

116.306 

147.711 

176.671 

120.639 

152.993 

182.712 

123.108 

155.991 

186.114 

123.637 

156.608 

186.756 

122.199 

154.798 

184.559 

0.93000 

48.428 

84.616 

118.814 

150.556 

179.487 

0.94000 

46.214 

80.808 

113.538 

143.912 

171.544 

0.95000 

43.298 

75.741 

106.451 

134.932 

160.772 

0.96000 

39.736 

69.500 

97.652 

123.708 

147.258 

0.97000 

35.575 

62.166 

87.244 

110.362 

131.132 

0.98000 

30.846 

53.804 

75.333 

95.034 

112.574 

0.99000 

25.576 

44.473 

62.024 

77.891 

91.808 

1.00000 

19.789 

34.234 

47.435 

59.120 

69.100 

SHEAR ( 7  ) 

SHEAR { 8 ) 

SHEAR ( 9 )  SHEAR ( 10 )  SUR.  ACCE 

-46.136 

-49.892 

-52.768 

-54.841 

0.55577 

-66.294 

-72.327 

-77.243 

-81.233 

0.59777 

-85.205 

-93.550 

-100.655 

-106.761 

0.61442 

102.696 

-113.435 

-122.861 

-131.092 

0.65476 

118.675 

-131.867 

-143.603 

-153.813 

0.71383 

133.071 

-148.657 

-162.524 

-174.535 

0.75987 

145.765 

-163.528 

-179.279 

-192.911 

0.77777 

156.573 

-176.196 

-193.578 

-208.624 

0.77939 

165.319 

-186.439 

-205.178 

-221.394 

0.78287 

171.929 

-194.128 

-213.889 

-231.001 

0.79422 

176.430 

-199.221 

-219.598 

-237.301 

0.80925 

178.881 

-201.756 

-222.304 

-240.253 

0.82491 

179.294 

-201  .819 

-222.107 

-239.948 

0.84298 

177.620 

-199.491 

-219.179 

-236.595 

0.86401 

173.800 

-194.828 

-213.706 

-230.467 

0.88189 

167.840 

-187.883 

-205.848 

-221.806 

0.88641 

159.842 

-178.753 

-195.723 

-210.761 

0.87013 

149.961 

-167.582 

-183.416 

-197.380 

0.83229 

138.341 

-154.512 

-168.999 

-181.676 

0.77709 

125.073 

-139.633 

-152.540 

-163.693 

0.70950 

110.199 

-122.977 

-134.105 

-143.541 

0.63296 

-93.767 

-104.571 

-113.774 

-121.371 

0.54988 

-75.873 

-84.495 

-91.655 

-97.332 

0.46294 

-56.675 

-62.914 

-67.902 

-71.575 

0.37529 

-36.390 

-40.057 

-42.719 

-44.287 

0.28939 

-15.265 

-16.190 

-16.373 

-15.736 

0.20557 

6.438 

8.398 

10.811 

13.709 

0.12150 

28.456 

33.385 

38.439 

43.594 

0.03298 

50.523 

58.414 

66.071 

73.416 

-0.06430 

72.364 

83.103 

93.243 

102.661 

-0.17299 

93.687 

107.067 

119.499 

130.825 

-0.29310 

114.185 

129.939 

144.416 

157.446 

-0.42223 

133.547 

151.372 

167.618 

182.116 

-0.55668 

151.465 

171.048 

188.772 

204.488 

-0.69270 

167.644 

188.673 

207.591 

224.274 

-0.82690 

181.805 

203.983 

223.822 

241.222 

-0.95590 

193.683 

216.735 

237.235 

255.109 

-1.07562 

203.033 

226.703 

247.617 

265.719 

-1.18115 

209.636 

233.670 

254.755 

272.844 

-1.26747 

213.301 

237.435 

258.444 

276.279 

-1.33067 

213.874 

237.817 

258.488 

275.826 

-1.36888 

211.242 

205.335 

196.127 

183.645 

167.968 

149.244 

127.684 

103.571 

77.244 


234.670 

227.892 

217.432 

203.304 

185.593 

164.466 

140.171 

113.040 

83.475 


254.718 

247.008 

235.292 

219.580 

199.971 

176.656 

149.919 

120.135 

87.756 


271.308 

262.592 

249.610 

232.382 

211.023 

185.753 

156.888 

124.837 

90.091 


-1.38236 

-1.37279 

-1.34222 

-1.29243 

-1.22447 

-1.13870 

-1.03484 

-0.91221 

-0.77001 


CYCLE  NUMBER  IS 


3 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  *  0 

TIME  SHEAR ( 2 )  SHEAR ( 3 )  SHEAR ( 4 )  SHEAR ( 5  ) 

SHEAR ( 6 ) 


1.01000 

12.665 

21.534 

29.222 

35.688 

40.877 

1.02000 

6.380 

10.565 

13.899 

16.252 

17.607 

1.03000 

-0.262 

-0.900 

-2.068 

-3.908 

-6.427 

1.04000 

-7.114 

-12.661 

-18.447 

-24.530 

-30.933 

1.05000 

-13.958 

-24.518 

-34.971 

-45.326 

-55.602 

1.06000 

-20.626 

-36.242 

-51.357 

-65.976 

-80.089 

1.07000 

-27.051 

-47.601 

-67.319 

-86.141 

-104.007 

1.08000 

-33.209 

-58.408 

-82.574 

-105.461 

-126.923 

1.09000 

-39.033 

-68.537 

-96.860 

-123.572 

-148.391 

1.10000 

-44.411 

-77.875 

-109.945 

-140.119 

-167.988 

1.11000 

-49.227 

-86.270 

-121.622 

-154.799 

-185.351 

1.12000 

-53.395 

-93.538 

-131.687 

-167.371 

-200.181 

1.13000 

-56.846 

-99.504 

-139.938 

-177.636 

-212.232 

1.14000 

-59.502 

-104.034 

-146.196 

-185.412 

-221.291 

1.15000 

-61.277 

-107.037 

-150.321 

-190.517 

-227.153 

1.16000 

-62.105 

-108.445 

-152.215 

-192.788 

-229.629 

1.17000 

-61.962 

-108.202 

-151.799 

-192.102 

-228.559 

1.18000 

-60.854 

-106.269 

-149.008 

-188.377 

-223.835 

1.19000 

-58.796 

-102.632 

-143.799 

-181.571 

-215.416 

1.20000 

-55.794 

-97.297 

-136.167 

-171.687 

-203.333 

1.21000 

-51 .850 

-90.291 

-126.157 

-158.787 

-187.691 

1.22000 

-46.976 

-81.658 

-113.861 

-143.005 

-168.668 

1.23000 

-41.198 

-71.470 

-99.417 

-124.548 

-146.519 

1.24000 

-34.569 

-59.837 

-83.011 

-103.676 

-121.559 

1.25000 

-27.176 

-46.920 

-64.875 

-80.695 

-94.155 

1.26000 

-19.136 

-32.921 

-45.285 

-55.945 

-64.709 

1.27000 

-10.598 

-18.081 

-24.557 

-29.801 

-33.659 

1.28000 

-1.729 

-2.667 

-3.033 

-2.672 

-1.478 

1.29000 

7.300 

13.045 

18.921 

24.998 

31.318 

1.30000 

16.325 

28.776 

40.927 

52.742 

64.183 

1.31000 

25.193 

44.253 

62.600 

80.073 

96.547 

1.32000 

33.763 

59.216 

83.558 

106.502 

127.832 

1.33000 

41.897 

73.412 

103.430 

131.552 

157.472 

1.34000 

49.463 

86.597 

121.864 

154.768 

184.923 

1.35000 

56.329 

98.538 

138.532 

175.729 

209.676 

1.36000 

62.367 

109.018 

153.128 

194.048 

231.267 

1.37000 

67.454 

117.833 

165.379 

209.382 

249.286 

1.38000 

71.484 

124.804 

175.039 

221.426 

263.374 

1.39000 

74.367 

129.778 

181.900 

229.926 

273.233 

1.40000 

76.032 

132.633 

185.795 

234.673 

278.622 

1.41000 

76.434 

133.283 

186.602 

235.517 

279.368 

1.42000 

75.544 

131.677 

184.247 

232.368 

275.372 

1.43000 

73.357 

127.805 

178.713 

225.203 

266.622 

1.44000 

69.886 

121.688 

170.035 

214.078 

253.192 

1.45000 

65.163 

113.392 

158.311 

199.122 

235.247 

1.46000 

59.246 

103.019 

143.695 

180.545 

213.045 

1.47000 

52.213 

90.716 

126.403 

158.626 

186.923 

1.48000 

44.172 

76.673 

106.707 

133.714 

157.295 

1.49000 

35.253 

61.121 

84.931 

106.215 

124.641 

1.50000 

25.612 

44.328 

61.442 

76.586 

89.499 

SHEAR ( 7 ) 

SHEAR { 8 ) 

SHEAR { 9 ) 

SHEAR(10 ) 

SUR.  ACCE 

44.763 

47.344 

48.621 

48.585 

-0.58820 

17.985 

17.414 

15.911 

13.475 

-0.39842 

-9.586 

-13.340 

-17.656 

-22.482 

-0.20729 

-37.622 

-44.543 

-51.618 

-58.707 

-0.02326 

-65.777 

-75.779 

-85.467 

-94.658 

0.16300 

-93.662 

-106.587 

-118.692 

-129.835 

0.35929 

120.832 

-136.481 

-150.817 

-163.752 

0.56127 

146.809 

-164.986 

-181.380 

-195.927 

0.75858 

171.118 

-191.641 

-209.914 

-225.876 

0.94538 

193.315 

-215.988 

-235.938 

-253.110 

1.12227 

212.996 

-237.573 

-258.972 

-277.134 

1.29047 

229.791 

-255.964 

-278.546 

-297.451 

1.44704 

243.377 

-270.773 

-294.217 

-313.586 

1.58591 

253.479 

-281.661 

-305.593 

-325.117 

1.70162 

259.868 

-288.354 

-312.353 

-331.701 

1.79120 

262.344 

-290.635 

-314.269 

-333.101 

1.85327 

260.741 

-288.354 

-311.206 

-329.196 

1.88679 

254.949 

-281.429 

-303.122 

-319.975 

1.89099 

244.933 

-269.857 

-290.057 

-305.516 

1.86558 

230.753 

-253.723 

-272.132 

-285.970 

1.81040 

212.558 

-233.205 

-249.542 

-261.553 

1.72469 

190.572 

-208.561 

-222.559 

-232.551 

1.60699 

165.090 

-180.116 

-191.531 

-199.327 

1.45619 

136.463 

-148.255 

-156.875 

-162.333 

1.27274 

105.099 

-113.423 

-119.080 

-122.097 

1.05915 

-71.459 

-76.128 

-78.702 

-79.218 

0.81959 

-36.046 

-36.938 

-36.355 

-34.349 

0.55926 

0.595 

3.536 

7.296 

11.809 

0.28375 

37.881 

44.648 

51.553 

58.518 

-0.00121 

75.200 

85.730 

95.696 

105.014 

-0.29017 

111.917 

126.093 

138.990 

150.525 

-0.57819 

147.386 

165.041 

180.699 

194.282 

-0.86091 

180.968 

201  .878 

220.088 

235.526 

-1.13420 

212.044 

235.928 

256.443 

273.515 

-1.39377 

240.031 

266.551 

289.081 

307.535 

-1.63491 

264.396 

293.155 

317.362 

336.920 

-1.85270 

284.666 

315.210 

340.713 

361.064 

-2.04234 

300.430 

332.258 

358.637 

379.451 

-2.19958 

311.348 

343.923 

370.732 

391.664 

-2.32090 

317.153 

349.920 

376.704 

397.406 

-2.40354 

317.662 

350.062 

376.366 

396.495 

-2.44541 

140 

312.779 

344.265 

369.650 

388.876 

-2.44505 

302.504 

332.549 

356.597 

374.608 

-2.40160 

286.938 

315.041 

337.359 

353.867 

-2.31482 

266.281 

291.972 

312.196 

326.935 

-2.18518 

240.830 

263.675 

281.467 

294.197 

-2.01392 

210.972 

230.574 

245.632 

256.141 

-1.80316 

177.176 

193.185 

205.244 

213.352 

-1.55603 

139.984 

152.104 

160.941 

166.504 

-1.27668 

100.003 

107.995 

113.439 

116.353 

-0.97027 

141 
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CYCLE  NUMBER  IS  4 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  =  0 


f 

9 


TIME  SHEAR { 2 )  SHEAR ( 3  )  SHE AR ( 4  )  SHEAR ( 5  ) 

SHE AR ( 6  ) 


1.51000 

1.52000 

1.53000 

1.54000 

1.55000 

1.56000 

1.57000 

1.58000 

1.59000 

1.60000 

1.61000 

1.62000 

1.63000 

1.64000 

1.65000 

1.66000 

1.67000 

1.68000 

1.69000 

1.70000 

1.71000 

1.72000 

1.73000 

1.74000 

1.75000 

1.76000 

1.77000 

1.78000 

1.79000 

1.80000 

1.81000 

1.82000 

1.83000 

1.84000 

1.85000 

1.86000 

1.87000 

1.88000 

1.89000 

1.90000 

1.91000 

1.92000 

1.93000 


13.830 

4.201 

-5.493 

-14.854 

-23.639 

-31.760 

-39.243 

-46.156 

-52.558 

-58.468 

-63.859 

-68.655 

-72.739 

-75.967 

-78.192 

-79.298 

-79.225 

-77.974 

-75.583 

-72.091 

-67.518 

-61.863 

-55.125 

-47.330 

-38.555 

-28.937 

-18.670 

-7.984 

2.879 

13.692 

24.265 

34.456 

44.165 

53.318 

61.849 

69.678 

76.702 

82.793 

87.812 

91.617 

94.088 

95.130 

94.681 


23.016 

7.078 

-9.077 

-25.048 

-40.449 

-54.979 

-68.482 

-80.933 

-92.379 

-102.849 

-112.299 

-120.593 

-127.531 

-132.905 

-136.540 

-138.319 

-138.168 

-136.043 

-131.918 

-125.780 

-117.642 

-107.549 

-95.577 

-81.841 

-66.514 

-49.831 

-32.096 

-13.659 

5.115 

23.878 

42.326 

60.208 

77.312 

93.458 

108.473 

122.180 

134.389 

144.896 

153.487 

159.955 

164.112 

165.805 

164.926 


30.956 

9.365 

-12.529 

-34.360 

-55.750 

-76.326 

-95.773 

-113.864 

-130.460 

-145.464 

-158.760 

-170.193 

-179.580 

-186.750 

-191.556 

-193.867 

-193.555 

-190.500 

-184.608 

-175.832 

-164.193 

-149.782 

-132.760 

-113.359 

-91.871 

-68.645 

-44.070 

-18.557 

7.474 

33.609 

59.454 

84.634 

108.804 

131.636 

152.824 

172.074 

189.103 

203.638 

215.415 

224.188 

229.731 

231.855 

230.417 


37.666 

10.994 

-16.047 

-43.112 

-69.852 

-95.901 

-120.869 

-144.353 

-165.956 

-185.334 

-202.233 

-216.496 

-228.029 

-236.744 

-242.507 

-245.142 

-244.459 

-240.292 

-232.525 

-221.112 

-206.091 

-187.599 

-165.880 

-141.274 

-114.185 

-85.048 

-54.311 

-22.420 

10.177 

43.020 

75.634 

107.530 

138.219 

167.225 

194.097 

218.419 

239.813 

257.936 

272.480 

283.171 

289.766 

292.059 

289.893 


43.158 

11.923 

-19.742 

-51.513 

-83.055 

-113.991 

-143.861 

-172.110 

-198.150 

-221.464 

-241.697 

-258.652 

-272.222 

-282.305 

-288.761 

-291.415 

-290.081 

-284.597 

-274.850 

-260.816 

-242.576 

-220.336 

-194.413 

-165.210 

-133.188 

-98.826 

-62.609 

-25.020 

13.448 

52.279 

90.916 

128.773 

165.242 

199.717 

231.621 

260.422 

285.643 

306.871 

323.749 

335.973 

343.287 

345.483 

342.410 


1.34000 

92.714 

161.418 

225.334 

283.172 

333.983 

1.95000 

89.230 

155.275 

216.586 

271.876 

320.204 

1.96000 

84.258 

146.539 

204.228 

256.074 

301.174 

1.97000 

77  .851 

135.301 

188.387 

235.935 

277.102 

1.98000 

70.084 

121.701 

169.268 

211.720 

248.302 

1.99000 

61.059 

105.924 

147.143 

183.784 

215.186 

2.00000 

50.905 

88.207 

122.353 

152.554 

178.248 

SHEAR ( 7  ) 

SHEAR { 8 ) 

SHEAR ( 9  )  SHE AR ( 10 )  SUR.  ACCE 

47.442 

50.522 

52.362 

52.898 

-0.58964 

12.118 

11.536 

10.115 

7.700 

-0.26461 

-23.709 

-28.043 

-32.834 

-38.077 

0.05406 

-59.723 

-67  .866 

-75.947 

-83.630 

0.34754 

-95.567 

-107.459 

-118.532 

-128.368 

0.61227 

-130.782 

-146.205 

-159.942 

-171.809 

0.85469 

-164.784 

-183.449 

-199.650 

-213  .459 

1  .08293 

-196.923 

-218.602 

-237.181 

-252.782 

1.30169 

-226.584 

-251.152 

-272.011 

-289.217 

1.51196 

-253.249 

-280.592 

-303.566 

-322.185 

1.71322 

-276.503 

-306.399 

-331.258 

-351.089 

1.90466 

-296.010 

-328.069 

-354.525 

-375.322 

2.08454 

-311.498 

-345.185 

-372.849 

-394.302 

2.24865 

-322.767 

-357.437 

-385.797 

-407.530 

2.38997 

-329.667 

-364.613 

-393.058 

-414.642 

2.50032 

-332.066 

-366.576 

-394.473 

-415.453 

2.57304 

-329.827 

-363.250 

-390.024 

-409.973 

2.60516 

-322.833 

-354.625 

-379.811 

-398.379 

2.59762 

-311.023 

-340.757 

-364.024 

-380.957 

2.55342 

-294.436 

-321.782 

-342.913 

-358.040 

2.47486 

-273.225 

-297.923 

-316.764 

-329.958 

2.36163 

-247.657 

-269.481 

-285.887 

-297.031 

2.21089 

-218.087 

-236.816 

-250.621 

-259.587 

2.01928 

-184.938 

-200.332 

-211.341 

-217.990 

1.78535 

-148.672 

-160.470 

-168.474 

-172.674 

1.51112 

-109.779 

-117.718 

-122.517 

-124.159 

1.20208 

-68.768 

-72.614 

-74.043 

-73.054 

0.86613 

-26.172 

-25.749 

-23.698 

-20.050 

0.51230 

17.446 

22.235 

27.806 

34.095 

0.14964 

61.488 

70.653 

79.712 

88.570 

-0.21368 

105.318 

118.790 

131.234 

142.538 

-0.57114 

148.267 

165.905 

181.572 

195.156 

-0.91848 

189.640 

211.244 

229.925 

245.589 

-1.25358 

228.740 

254.050 

275.501 

293.020 

-1.57559 

264.894 

293.588 

317.530 

336.659 

-1.88351 

297.471 

329.160 

355.273 

375.745 

-2.17499 

325.907 

360.129 

388.044 

409.563 

-2.44580 

349.712 

385.934 

415.222 

437.458 

-2.69009 

368.473 

406.098 

436.278 

458.870 

-2.90124 

381.848 

420.238 

450.787 

473.350 

-3.07286 

389.565 

428.062 

458.443 

480.587 

-3.19946 

391.421 

429.377 

459.065 

480.412 

-3.27684 

T I  ME  = 


2.0000  SEC 


CYCLE  NUMBER  IS 


5 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  =  0 

TIME  SHEAR { 2  )  SHEAR ( 3  )  SHEAR ( 4  )  SHEAR ( 5 


SHEAR  ( 6 ) 

2.01000  34 
2.02000  23 
2.03000  11 
2.04000  -0 
2.05000  -11 
2.06000  -20 
2.07000  -28 
2.08000  -36 
2.09000  -43 
2.10000  -49 
2.11000  -56 
2.12000  -62 
2.13000  -69 
2.14000  -75 
2.15000  -82 
2.16000  -87 
2.17000  -91 
2.18000  -94 
2.19000  -96 
2.20000  -96 
2.21000  -95 
2.22000  -93 
2.23000  -90 
2.24000  -87 
2.25000  -82 
2.26000  -76 
2.27000  -68 
2.28000  -59 
2.29000  -49 
2.30000  -38 
2.31000  -25 
2.32000  -12 
2.33000  0 
2.34000  13 
2.35000  26 
2.36000  38 
2.37000  50 
2.38000  61 
2.39000  72 
2.40000  82 
2.41000  91 
2.42000  100 


.390  55.712 
.266  38.549 
.373  20.413 
.374  1.927 
.168  -16.193 
.621  -33.273 
.836  -48.882 
.203  -62.985 
.109  -75.920 
.793  -88.185 
.373  -100.183 
.927  -112.051 
.490  -123.644 
.978  -134.623 
.133  -144.578 
.568  -153.117 
.894  -159.929 
.850  -164.818 
.368  -167.696 
.552  -168.557 
.593  -167.433 
.660  -164.351 
.834  -159.311 
.072  -152.274 
.237  -143.165 
.156  -131.888 
.689  -118.366 
.782  -102.594 
.498  -84.694 
.018  -64.937 
.615  -43.732 
.624  -21.571 
.600  1.040 
.734  23.645 
.521  45.877 
.799  67.467 
.494  88.235 
.599  108.060 
.132  126.846 
.090  144.484 
.420  160.830 
.008  175.685 


73.361  87 

51.236  61 

28.162  34 

4.519  6 

-19.257  -21 

-42.683  -50 

-65.289  -79 

-86.734  -107 

-106.898  -135 

-125.872  -161 

-143.827  -185 

-160.830  -207 

-176.745  -227 

-191.273  -244 

-204.090  -259 

-214.956  -272 

-223.729  -282 

-230.295  -290 

-234.494  -295 

-236.108  -297 

-234.906  -296 

-230.704  -290 

-223.392  -280 

-212.931  -266 

-199.329  -248 

-182.636  -226 

-162.957  -201 

-140.474  -172 

-115.457  -141 

-88.270  -107 

-59.351  -71 

-29.187  -34 

1.725  3 

32.906  42 

63.920  80 

94.376  119 

123.931  157 

152,269  193 

179.098  227 

204.132  259 

227.091  288 

247.689  314 


.997  100.250 
.641  70.116 
.408  39.134 
.495  7.333 
.910  -25.302 
.596  -58.764 
.297  -92.887 
.644  -127.207 
.157  -160.921 
.285  -193.039 
.519  -222.663 
.510  -249.237 
.129  -272.617 
.424  -292.944 
.492  -310.420 
.354  -325.126 
.889  -336.961 
.851  -345.667 
.920  -350.889 
.753  -352.230 
.014  -349.296 
.406  -341.745 
.707  -329.348 
.819  -312.043 
.800  -289.965 
.865  -263.440 
.353  -232.932 
.675  -198.969 
.267  -162.071 
.558  -122.698 
.960  -81.246 
.876  -38.076 
.280  6.437 
.063  51.855 
.994  97.652 
.554  143.210 
.210  187.849 
.425  230.857 
.686  271.537 
.517  309.238 
.482  343.382 
.196  373.479 
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2.43000 

107.682 

188.803 

265.630 

336.313 

399.117 

2.44000 

114.239 

199.901 

280.612 

354.523 

419.954 

2.45000 

119.472 

208.693 

292.333 

368.538 

435.696 

2.46000 

123.198 

214.910 

300.507 

378.087 

446.072 

2.47000 

125.277 

218.335 

304.877 

382.908 

450.827 

2.48000 

125.619 

218.807 

305.229 

382.760 

449.719 

2.49000 

124.184 

216.235 

301.407 

377.439 

442.546 

2.50000 

120.971 

210.581 

293.314 

366.800 

429.166 

SHEAR ( 7 ) 

SHEAR { 8 ) 

SHEAR { 9 ) 

SHE AR ( 10  ) 

SUR.  ACCE 

110.598 

119.281 

126.232 

130.910 

-1.50824 

76.903 

82.032 

85.164 

85.310 

-1.00188 

42.267 

43.523 

42.258 

37.958 

-0.52648 

6.554 

3.538 

-2.213 

-9.691 

-0.14898 

-30.379 

-37.856 

-47.394 

-56.882 

0.14975 

-68.491 

-80.133 

-92.360 

-103.395 

0.41624 

-107.372 

-122.462 

-136.566 

-149.146 

0.67391 

-146.248 

-164.060 

-179.816 

-193.965 

0.91960 

-184.167 

-204.393 

-221.973 

-237.579 

1.14824 

-220.258 

-243.100 

-262.731 

-279.651 

1.36921 

-253.865 

-279.779 

-301.599 

-319.771 

1.60072 

-284.530 

-313.877 

-337.980 

-357.403 

1.85526 

-311.916 

-344.748 

-371.209 

-391.860 

2.13227 

-335.775 

-371.799 

-400.568 

-422.333 

2.42080 

-355.962 

-394.573 

-425.332 

-447.959 

2.70495 

-372.420 

-412.757 

-444.850 

-467.921 

2.96690 

-385.103 

-426.138 

-458.638 

-481.570 

3.18880 

-393.876 

-434.548 

-466.421 

-488.524 

3.35690 

-398.468 

-437.831 

-468.123 

-488.725 

3.46684 

-398.512 

-435.835 

-463.827 

-482.405 

3.52545 

-393.643 

-428.430 

-453.723 

-469.990 

3.54625 

-383.591 

-415.560 

-438.056 

-451.955 

3.54088 

-368.244 

-397.280 

-417.089 

-428.711 

3.51220 

-347.664 

-373.763 

-391.066 

-400.544 

3.45298 

-322.074 

-345.270 

-360.204 

-367.623 

3.35009 

-291.832 

-312.112 

-324.696 

-330.060 

3.19064 

-257.379 

-274.617 

-284.738 

-287.977 

2.96645 

-219.198 

-233.129 

-240.568 

-241.572 

2.67554 

-177.768 

-188.025 

-192.500 

-191.157 

2.32160 

-133.539 

-139.735 

-140.960 

-137.175 

1.91270 

-86.935 

-88.743 

-86.486 

-80.205 

1.46040 

-38.376 

-35.585 

-29.718 

-20.943 

0.97871 

11.694 

19.155 

28.633 

39.819 

0.48254 

62.771 

74.852 

87.812 

101.241 

-0.01454 

114.264 

130.831 

147.042 

162.472 

-0.50246 

165.497 

186.362 

205.542 

222.687 

-0.97609 

215.712 

240.666 

262.527 

281.098 

-1.43494 

264.107 

292.926 

317.210 

336.947 

-1.88142 

309.877 

342.318 

368.805 

389.481 

-2.31827 

352.257 

388.038 

416.521 

437.935 

-2.74611 

390.558 

429.334 

459.578 

481.523 

-3.16203 

424.187 

452.644 

475.514 

492.445 

503.139 

507.337 

504.830 

495.469 

479.183 


465.518 

495.984 

520.208 

537.749 

548.258 

551.477 

547.248 

535.510 

516.306 


497.221 

528.742 

553.511 

571.005 

580.831 

582.741 

576.631 

562.537 

540.611 


519.456 

550.971 

575.384 

592.129 

600.801 

601.167 

593.173 

576.919 

552.638 


-3.55934 

-3.92856 

-4.25892 

-4.53996 

-4.76275 

-4.92063 

-5.00930 

-5.02655 

-4.97173 


T I  ME  = 


2.5000  SEC 


INITIAL  EFFECTIVE  STRESS  NEW  PORE  PRESSURE 


239.3030 

0.0000 

627.6061 

378.4290 

925.6061 

668.7290 

1223.6061 

949.3483 

1521.6061 

1212.9518 

1819.6061 

1448.7770 

2117.6061 

1644.7080 

2415.6061 

1791.8976 

2713.6061 

1887.0133 

3011.6061 

1931.4199 

CYCLE  NUMBER  IS 


6 


DAMPING  RATIO  =  0.02000 

LAYER  FIRST  TO  LIQUEFY  =  0 

TIME  SHEAR ( 2 )  SHEAR ( 3 )  SHEAR ( 4 )  SHEAR ( 5 ) 

SHEAR ( 6 ) 


2.51000 

87.577 

123.660 

140.829 

150.300 

161.236 

2.52000 

76.665 

113.336 

132.704 

143.587 

154.610 

2.53000 

60.781 

98.666 

122.468 

136.672 

148.755 

2.54000 

43.384 

80.875 

109.946 

129.030 

143.165 

2.55000 

27.904 

61.835 

95.253 

120.078 

137.158 

2.56000 

16.417 

43.790 

78.985 

109.319 

129.877 

2.57000 

9.183 

28.739 

62.229 

96.476 

120.344 

2.58000 

5.240 

17.814 

46.339 

81.601 

107.592 

2.59000 

3.348 

11.010 

32.517 

65.068 

90.890 

2.60000 

2.592 

7.399 

21.392 

47.467 

69.987 

2.61000 

2.439 

5.618 

12.784 

29.416 

45.296 

2.62000 

2.496 

4.342 

5.776 

11.379 

17.900 

2.63000 

2.327 

2.510 

-0.968 

-6.434 

-10.641 

2.64000 

1.432 

-0.657 

-8.771 

-24.048 

-38.655 

2.65000 

-0.654 

-5.729 

-18.590 

-41.605 

-64.753 

2.66000 

-4.278 

-13.107 

-30.837 

-59.228 

-88.129 

2.67000 

-9.619 

-22.970 

-45.393 

-76.925 

-108.660 

2.68000 

-16.675 

-35.198 

-61.750 

-94.580 

-126.801 

2.69000 

-25.233 

-49.336 

-79.192 

-112.018 

-143.331 

2.70000 

-34.859 

-64.630 

-96.962 

-129.080 

-159.075 

2.71000 

-44.922 

-80.146 

-114.371 

-145.685 

-174.689 

2.72000 

-54.675 

-94.932 

-130.863 

-161.821 

-190.561 

2.73000 

-63.391 

-108.169 

-146.032 

-177.504 

-206.814 

2.74000 

-70.501 

-119.284 

-159.611 

-192.711 

-223.367 

2.75000 

-75.698 

-128.007 

-171.462 

-207.334 

-239.997 

2.76000 

-78.983 

-134.369 

-181.546 

-221.170 

-256.378 

2.77000 

-80.636 

-138.657 

-189.905 

-233.942 

-272.086 

2.78000 

-81.124 

-141.329 

-196.647 

-245.338 

-286.601 

2.79000 

-80.990 

-142.918 

-201.927 

-255.052 

-299.322 

2.80000 

-80.729 

-143.931 

-205.922 

-262.809 

-309.610 

2.81000 

-80.697 

-144.769 

-208.800 

-268.377 

-316.862 

2.82000 

-81.067 

-145.658 

-210.675 

-271.571 

-320.594 

2.83000 

-81.820 

-146.625 

-211.571 

-272.254 

-320.517 

2.84000 

-82.778 

-147.499 

-211.392 

-270.338 

-316.572 

2.85000 

-83.650 

-147.941 

-209.923 

-265.790 

-308.930 

2.86000 

-84.088 

-147.505 

-206.857 

-258.637 

-297.956 

2.87000 

-83.731 

-145.704 

-201.853 

-248.961 

-284.135 

2.88000 

-82.254 

-142.089 

-194.602 

-236.891 

-267.992 

2.89000 

-79.400 

-136.311 

-184.889 

-222.598 

-250.012 

2.90000 

-75.014 

-128.178 

-172.642 

-206.275 

-230.582 

2.91000 

-69.059 

-117.685 

-157.950 

-188.122 

-209.958 

2.92000 

-61.633 

-105.019 

-141.057 

-168.344 

-188.260 

2.93000 

-52.959 

-90.533  - 

122.323  - 

147.135  - 

165.497 

2.94000 

-43.366 

-74.709  - 

102.189  - 

124.682  - 

141.614 

2.95000 

-33.248 

-58.096 

-81.121  - 

101.168  - 

116.535 

2.96000 

-23.017 

-41.241 

-59.571 

-76.780 

-90.218 

2.97000 

-13.044 

-24.634 

-37.941 

-51.711 

-62.691 

2.98000 

-3.619 

-8.659 

-16.563 

-26.176 

-34.081 

2.99000 

5.082 

6.440 

4.315 

-0.408 

-4.613 

3.00000 

13.006 

20.562 

24.523 

25.335 

25.387 

SHEAR ( 7 ) 

SHEAR ( 8 ) 

SHEAR ( 9 )  SHEAR ( 10 )  SUR.  ACCE 

178.407 

201.927 

230.108 

259.839 

-4.26727 

170.561 

191.857 

216.031 

235.155 

-3.13912 

163.643 

182.031 

199.555 

204.937 

-2.14779 

157.125 

171.370 

179.361 

172.945 

-1.53820 

150.167 

158.622 

155.633 

141.880 

-1.20646 

141.635 

142.845 

129.842 

112.990 

-0.99159 

130.287 

123.794 

103.850 

86.413 

-0.83713 

115.072 

101.978 

79.005 

61.769 

-0.76275 

95.440 

78.335 

55.733 

38.619 

-0.77604 

71.520 

53.758 

33.662 

16.615 

-0.83998 

44.093 

28.742 

12.038 

-4.570 

-0.90217 

14.385 

3.313 

-9.887 

-25.368 

-0.92991 

-16.224 

-22.764 

-32.662 

-46.359 

-0.91405 

-46.462 

-49.678 

-56.620 

-68.209 

-0.85035 

-75.329 

-77.340 

-81.933 

-91.539 

-0.72551 

-102.166 

-105.327 

-108.629 

-116.777 

-0.52027 

-126.649 

-133.004 

-136.559 

-144.047 

-0.22164 

-148.762 

-159.722 

-165.344 

-173.103 

0.16851 

-168.759 

-184.999 

-194.372 

-203.313 

0.63325 

-187.108 

-208.607 

-222.865 

-233.701 

1.14515 

-204.397 

-230.561 

-249.990 

-263.056 

1.67082 

-221.221 

-251.024 

-274.985 

-290.093 

2.17454 

-238.047 

-270.192 

-297.245 

-313.644 

2.62279 

-255.107 

-288.178 

-316.369 

-332.849 

2.99028 

-272.335 

-304.948 

-332.154 

-347.280 

3.26562 

-289.377 

-320.293 

-344.576 

-356.979 

3.45360 

-305.646 

-333.868 

-353.752 

-362.402 

3.57272 

-320.420 

-345.254 

-359.907 

-364.282 

3.64901 

-332.948 

-354.041 

-363.331 

-363.458 

3.70869 

-342.538 

-359.901 

-364.335 

-360.712 

3.77218 

-348.642 

-362.642 

-363.206 

-356.641 

3.85070 

-350.899 

-362.211 

-360.159 

-351.591 

3.94564 

-349.164 

-358.673 

-355.313 

-345.649 

4.04993 

-343.512 

-352.171 

-348.670 

-338.676 

4.15045 

-334.206 

-342.877 

-340.132 

-330.375 

4.23078 

-321.652 

-330.959 

-329.526 

-320.362 

4.27367 

-306.335 

-316.566 

-316.646 

-308.244 

4.26298 

-288.749 

-299.832 

-301.299 

-293.673 

4.18523 

-269.333 

-280.884 

-283.355 

-276.402 

4.03078 

-248.429 

-259.853 

-262.771 

-256.311 

3.79498 

-226.257 

-236.881 

-239.618 

-233.432 

3.47889 

202.918 

178.419 

152.706 

125.711 

-97.393 

-67.778 

-36.981 

-5.223 

27.174 


-212.116 

-185.714 

-157.827 

-128.608 

-98.215 

-66.825 

-34.640 

-1.903 

31.102 


-214.076 

-186.421 

-157.003 

-126.208 

-94.433 

-62.056 

-29.419 

3.185 

35.516 


-207.948 

-180.183 

-150.573 

-119.627 

-87.870 

-55.802 

-23.846 

7.677 

38.558 


3.08973 

2.64050 

2.14871 

1.63431 

1.11715 

0.61438 

0.13842 

-0.30408 

-0.71230 
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239.3030 

627.6061 

925.6061 

1223.6061 

1521.6061 

1819.6061 

2117.6061 

2415.6061 

2713.6061 

3011.6061 


0.0000 

444.3646 

808.9455 

1210.6247 

1634.8356 

2022.8055 

2305.7583 

2454.2316 

2493.3637 

2459.3924 


END 

FILMED 

3-85 


DTIC 


